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Abstract. An analysis is presented of VLT-FLAMES spectroscopy for three Galactic clusters, NGC 3293, NGC 4755 and 
NGC 661 1. Non-LTE model atmosphere calculations have been used to estimate effective temperatures (from either the helium 
spectrum or the silicon ionization equilibrium) and gravities (from the hydrogen spectrum). Projected rotational velocities have 
been deduced from the helium spectrum (for fast and moderate rotators) or the metal line spectrum (for slow rotators). The 
origin of the low gravity estimates for apparently near main sequence objects is discussed and is related to the stellar rotational 
velocity. The atmospheric parameters have been used to estimate cluster distances (which are generally in good agreement 
with previous determinations) and these have been used to estimate stellar luminosities and evolutionary masses. The observed 
Hertzsprung-Russell diagrams are compared with theoretical predictions and some discrepancies including differences in the 
main sequence luminosities are discussed. Cluster ages have been deduced and evidence for non-coeval star formation is found 
for all three of the clusters. Projected rotational velocities for targets in the older clusters, NGC 3293 and NGC 4755, have 
been found to be systematically larger than those for the field, confirming recent results in other similar age clusters. The 
distribution of projected rotational velocities are consistent with a Gaussian distribution of intrinsic rotational velocities. For 
the relatively unevolved targets in the older clusters, NGC 3293 and NGC 4755, the peak of the velocity distribution would be 
250 km s~^with a full-width-half-maximum of approximately 180 km s~^. For NGC 6611, the sample size is relatively small 
but implies a lower mean rotational velocity. This may be evidence for the spin-down effect due to angular momentum loss 
through stellar winds, although our results are consistent with those found for very young high mass stars. For all three clusters 
we deduce present day mass functions with F-values in the range of -1.5 to -1.8, which are similar to other young stellar clusters 
in the Milky Way. 

Key words, stars: early-type - stars: fundamental parameters - stars: rotation - Hertzsprung-Russell (HR) and C-M diagrams 
- open clusters and associations: NGC 3293, NGC 4755 & NGC 661 1 



1. Introduction 

As part of a European Southern Observatory Large 
Programme, we are using the Fibre Large Array Multi- 
Element Spectrograph (FLAMES) at the Very Large Telescope 
to survey the hot stellar populations in our Galaxy and in 
the Magellanic Clouds. We have observed in excess of 50 
O-type stars and 500 B-type stars, in a total of seven clusters. 
These data have been supplemented with observations of 
our brighter targets using the Fibre-Fed Extended Range 
Optical Spectrograph (FEROS) at La Silla. Evans et al. ("2005 
hereafter Paper I) have given an overview of the scientific 
goals of the project, discussed the target selection and data 
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reduction techniques, and presented some preliminary analysis 
for three Galactic clusters. The scientific aims of the survey 
as discussed in Paper I, included understanding how rotation 
and metallicity influence stellar evolution and stellar winds, 
calibration of the wind momentum-luminosity relationship 
and the nature of supernova progenitors. In the case of stellar 
evolution, theoretical studies (see, for example, Heger & 
Langer l20001 Maeder & Meynet ;200()| provide predictions of 
the surface enhancements expected for helium and nitrogen. 
Of particular relevance to our survey are the larger logarithmic 
enhancements predicted for lower metallicity regimes (Maeder 
& Meynet '2001 ). Although there are many previous detailed 
quantitative studies of O- and B-type stars (Crowther at al. 
I2OO2I I2OO6I Bouret et al. l2003l Hillier at al. l2003l Evans et 
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al. l2004l Tr undle et al. l2004l W alborn et al. l2004l Dufton 
et al. 120051 Trundle & Lennon I2005> . these have tended 
to concentrate on narrow lined (and hence possibly slowly 
rotating) stars, whilst the sample sizes have been relatively 
small. The current survey will sample both the main sequence 
and supergiant regimes and include rapidly rotating stars, 
thereby extending the scope of such studies. 

In order to have a sample of stars in the Milky Way to 
compare with the FLAMES datasets for the SMC and LMC 
we have obtained spectra of likely cluster members in the 
three Galactic Clusters, NGC 3293, NGC 4755 and NGC 66 11 . 
These clusters were chosen so as to have a similar age 
range to their lower metallicity Magellanic Cloud counter- 
parts. Although the Milky Way has an abundance gradient 
(see, for example, Afflerbach et al. [19971 Rolleston et al. 2000, 
Gummersbach et al. 1998, Daflon & Cunha 2004) these clus- 
ters are relatively close to the solar position and appear to have 
a similar metallicity to that of the Sun from a preliminary model 
atmosphere analysis of selected narrow line targets (see Hunter 
et al. 2006 for details). Previous studies have discussed the 
HR-diagrams, star formation timescales, initial mass funtions 
and rotational velocity distributions in young Galactic clus- 
ters and have highlighted the value of spectroscopic analysis 
of cluster members rather than spectral type assignments from 
photometric colours (e.g. see Massey et al. 119951 Hillenbrand 
et al. I1993> . Additionally there has been an revival of inter- 
est in the rotational velocity distribution of large samples of 
massive stars for comparison with the predictions of theoret- 
ical evolutionary models and also to probe initial conditions 
during star-formation. Strom et al. C2005'\ have measured pro- 
jected rotational velocities for approximately 200 B-type stars 
in h&xPer and suggest that differences between the mean val- 
ues for the cluster and field stars are due to the star-formation 
process rather than being a consequence of angular momen- 
tum evolution as they have moved from the ZAMS. Huang & 
Gies (1200 5^ have presented rotational velocities of stars in nine 
Galactic open clusters with age ranges of approximately 6 to 73 
Myr, and have provided an interpretation on the basis of evolu- 
tionary arguments. Additionally the extensive field star sample 
of Abt et al. (2QQ2 1 remains a benchmark study with which to 
compare cluster rotational velocities. All of these studies illus- 
trate the importance of large stellar samples in order to amelio- 
rate statistical significance problems arising from the random 
incUnation angle of the rotational axes. The availability of the 
FLAMES spectrograph on an 8m telescope has allowed stud- 
ies of large samples to be extended to the Magellanic Clouds. 
For example, Martayan et al. J2006t have studied the effects of 
metallicity on the rotational velocities of B-type and Be stars in 
the LMC and find that the LMC stars rotate faster than Be stars 
at solar metallicity. 

In this paper we analyse the early type stars in the three 
Galactic clusters observed with the FLAMES spectrograph. 
Stellar atmospheric parameters are deduced using non-LTE 
model atmosphere calculations and rotational velocities are 
estimated by comparing observed and theoretical rotationally 
broadened spectra. The former are discussed with particular 
reference to the surface gravities found for near main sequence 
objects; evidence is presented that relates the lower gravity es- 



timates to large stellar rotational velocities. Our projected rota- 
tional velocities are discussed both in the context of possible 
systematic differences between field and cluster populations 
and to infer the underlying rotational velocity distributions. The 
estimates of the atmospheric parameters are used to deduce 
cluster distances leading to estimates of stellar luminosities 
and evolutionary masses. The former are then used to generate 
Hertzsprung-Russell diagrams and the agreement with theoret- 
ical predictions is discussed, while present day mass functions 
are estimated from the latter 

2. Observational data 

The observational data have been obtained during a European 
Southern Observatory (ESO) Large Programme to study early- 
type stars in our Galaxy and the Magellanic Clouds. For the 
former, targets were observed in three clusters - NGC 3293, 
NGC 4755 and NGC 661 1 - with the Fibre Large AiTay Multi- 
Element Spectrograph (FLAMES, R « 20000, see Pasquini 
et al. 2002 1 on the Very Large Telescope (VLT). The targets 
were selected based on astrometry and photometry from the 
pre-FLAMES observations using the wide-field imager (WFI) 
on the 2.2-m Max Planck Gesellschaft (MPG)/ESO telescope, 
whilst the brightest stars were observed separately with the 
Fibre-Fed, Extended Range Optical Spectrograph (FEROS, R 
w 48 000, see Kaufer et al. .1999) on the MPG. Further details 
of the target selection criteria and data acquisition are presented 
in Paper I. Using six standard high-resolution GIRAFFE set- 
tings, the FLAMES data cover a wavelength range of 3950- 
4755A and 6380-6620A with a typical signal-to-noise (S/N) 
ratio of 100, whilst the FEROS spectra range from 3600 to 
9200 A with a S/N ratio >100. Reductions were performed 
using both IRAF' software and the FLAMES pipeline (gir- 
BLDRS, Blecha et al. 20031, which showed excellent agree- 
ment. Details of the data reduction procedure, spectral classi- 
fication and radial velocity measurements of the targets were 
presented in Paper I, which also contains the coordinates and 
cross-referencing for the designations of the targets. As an ap- 
pendix we include the finding charts for all objects we observed 
in these clusters (available on-line only). 

In this paper, we have only considered targets classified as 
earlier than B9. This cutoff was adopted as for later spectral 
types the neutral helium lines, which have been used as a tem- 
perature diagnostic were not normally visible. One exception 
is the target 4755-001 (note we have adopted the nomenclature 
used in Paper I), which despite a classification of B9 la had a 
measurable helium spectrum. 

3. Method of analysis 

In this section, we discuss the methods that were used to esti- 
mate the stellar atmospheric parameters, projected rotational 
velocities, masses and luminosities together with the cluster 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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distances. We also discuss some problematic stars (Sect. 13 .4> 
and also the criteria to identify cluster membership (Sect. ll!6t . 

3. 1. Initial estimation of atmospheric parameters 

The estimation of the stellar atmospheric parameters - effec- 
tive temperature (T^ff) and logarithmic surface gravity (log g) 
- requires an iterative process and we have adopted different 
approaches depending on the stellar effective temperature. For 
B2-B8 stars, which make up the majority of the sample, the 
strengths of the hydrogen and neutral helium lines were uti- 
Uzed, while the ionization balance of Si Ill/iv was used for 
the earlier B spectral types. In the case of the O-type stars, the 
hydrogen and neutral and ionized helium line profiles provided 
diagnostics for the atmospheric parameters, together with those 
of the wind. The different methodologies are discussed below: 

3.1 .1 . B2 - B8 Stars: H & He lines 

As the FLAMES dataset contains a relatively large number of 
mid and late B-type stars, we have utilized relatively simple 
methods to estimate their atmospheric parameters. The ratio 
of the strengths of the neutral helium and hydrogen lines is 
sensitive to effective temperature for spectral types later than 
Bl, whilst the strength of the hydrogen lines can be used to 
estimate the surface gravity (see for example Dufton et al. 
|1999) . For hydrogen, the equivalent widths of the H 7 and H 6 
lines were measured with the continuum defined at ± 16A from 
the line centre. For any given star these normally agreed to 
within 3% (which is consistent with our theoretical calcula- 
tions) and hence their average was adopted. For helium, while a 
number of different features were available, the best observed 
line which was seen in effectively all our spectra was that at 
4026A. Moreover, this line has been found to yield reliable ef- 
fective temperature estimates in a study of hot stars towards the 
Galactic centre by Dufton et al. ( 1999 1. The equivalent widths 
(EWs) of this line were estimated with the continuum normally 
defined at ±6A from the line centre, although for some of the 
sharper and weaker He I lines a narrower range was adopted 
where this was believed to increase the reliability of the mea- 
surements. 

Diagnostic plots (see Fig.[0 were deduced from synthetic 
spectra that had been calculated using the non-LTE model at- 
mosphere codes TLUSTY/SYNSPEC (Hubeny 1988, Hubeny & 
Lanz 1995i Hubeny et al. .l998 see also Ryans et al. 2003 ). 
A microturbulent velocity of 5 kms^^ was adopted, consis- 
tent with most of our targets having spectral types that imply 
that they are on and/or near the main sequence (see, for exam- 
ple, Gies & Lambert 1 1992> . However, as discussed by Dufton 
et al. J 1999). this choice is not crucial for these broad hydro- 
gen and diffuse helium lines. Fig. Q shows the variation of 
the relative strength of the neutral helium and hydrogen lines 
(EW(He/H)), and of the mean of the equivalent width of the hy- 
drogen lines (EW(H)), for the range of effective temperatures 
that were available from our grid of model atmospheres. Note 
that all measurements of the theoretical equivalent widths have 
used the same continuum definitions as for the observations. 



For estimating the gravity, the hydrogen lines provide reliable 
diagnostics across the entire temperature range. By contrast the 
helium line diagnostic flattens out at temperatures >20000 K, 
which sets the effective upper limit for which this method can 
be used. 

For each star, estimates of the atmospheric parameters were 
obtained by an iterative procedure, using the diagnostic plots 
illustrated in Fig. Firstly the observed value EW(He/H) 
was compared to theoretical predictions for a surface gravity 
of log g = 4.0 dex (g in cms~^) to give an initial estimate 
of the effective temperature (Fig. [^). The observed value of 
log(EW(H)) was then compared with theoretical predictions to 
give an estimate of surface gravity (Fig.fQ)). The above proce- 
dure was then repeated until convergence was found for both 
effective temperature and surface gravity. 

3.1 .2. B0.5 - Bl .5: Si ionization balance 

The neutral helium spectrum is not a useful effective temper- 
ature diagnosis for stars earlier than approximately B2 (Teg 
>20000 K, see Fig. [0 and hence for these objects, the sil- 
icon ionization balance was utilized instead. However, atmo- 
spheric parameters could not be estimated for stars with large 
rotational velocities where the silicon spectra was poorly ob- 
served. Therefore we only considered narrow-lined stars and 
adopted their temperature estimates as representative for other 
stars with the same spectral type. 

By adopting an appropriate surface gravity (4.0 dex for 
dwarfs and 3.0 dex for giants) and a microturbulence (5 
km s^^but see Sect. l3.4l for a discussion of the microturbulence 
of the supergiant targets) an initial effective temperature could 
be estimated using the relative strengths of the Si Ill/iv lines 
and the TLUSTY (Hubeny 1988 Hubeny & Lanz 1995 1 non- 
LTE model atmosphere grids discussed in Sect. B.l.Tl fsee also 
Dufton et al. 120051 for details on the calculation of the silicon 
spectrum). These temperature estimates were in turn used to es- 
timate surface gravity by the comparison of observed and theo- 
retical hydrogen Balmer line profiles. The new gravity estimate 
was then used as the starting point in the next iteration to es- 
timate the effective temperature and the process was repeated 
until convergence was obtained. 

Narrow-lined stars in the clusters were selected with spec- 
tral types between BO. 5 and B1.5 for both dwarfs and giants. 
For BO. 5 V, the target 6611-012 was analyzed, whilst for Bl 
V, four narrow-lined stars were available (3293-018,4755-015, 
6611-021 and 6611-033). For B1.5 V, the spectra of all three 
targets (4755-011, 4755-017 and 6611-032) only show Si III 
lines and the effective temperature for this spectral type was 
adopted as the average of those of Bl and B2 type stars. For 
Bl 111, the effective temperature estimates of three stars (3293- 
003, 3293-007 and 3293-010) were averaged. No narrow fined 
star was present in our sample with a spectral type B 1 .5 111. As 
in the case of dwarfs, for this spectral type, the effective tem- 
perature was taken to be the average of those for our B 1 111 and 
B2 111 stars. In the case of the latter spectral type, the targets 
3293-034 and 4755-003 were excluded from the estimation of 
mean effective temperature. This was because the former has a 



4 



Dufton et al: Stellar parameters in NGC 3293, NGC 4755 and NGC 661 1 



Fig. 1. Diagnostic plots used for the initial estimation of atmospheric parameters assuming a microturbulence of 5 km s~^: (left) 
the ratio of the equivalent width of the He I line at 4026A to the mean of the Hydrogen H 7 and H S equivalent widths; the ratio 
is plotted as a function of effective temperature (Tcs) for three logarithmic surface gravities:(rig/zf) the mean equivalent widths 
of the H 7 and H5 lines on a logarithmic scale plotted as a function of effective temperature (Toff) for five logarithmic surface 
gravities. 




Table 1. Methods used to estimate the effective temperature scale for the BO. 5-1. 5 spectral-types 



Spec. Type 




Method 


B0.5 V 


27 000 K 


Si III/IV ionization balance (6611-012) 


Bl.O V 


25 000 K 


Si III/IV ionization balance (3293-018; 4755-015; 6611-021 & -033) 


B1.5 V 


22 500 K 


Interpolation between B 1 V and B2 V 


B2.0 V 


20 000 K 


Average of the effective temperature estimate for 1 1 stars based on their helium spectra 


Bl.O III 


21 500 K 


Si III/IV ionization balance (3293-003, -007 & -010) 


B1.5 III 


20 500 K 


Interpolation between B 1 III and B2 III 


B2.0 III 


19 500 K 


Average of the effective temperature estimate for 2 stars based on their helium spectra 



pecuUar spectrum (see paper I), whilst the estimate for the lat- 
ter (Tcff=24000 K) was beyond the high temperature limit for 
using the helium spectra as a reliable diagnostic. 

These effective temperature estimates are summarized in 
Table ^ where they have been rounded to the nearest 500K. 
They were adopted for all other stars with the same spectral 
classification but with individual surface gravities being de- 
duced from the hydrogen Balmer line profiles. 

3.1.3. O and BO 

The BO and O-type stars have temperatures at the edge of, or 
outside, the range of our TLUSTY grid. In addition many of 
these objects show wind features in their optical spectra. For 
these reasons the hottest objects in our sample have been an- 
alyzed with the unified model atmosphere code FASTWIND as 
described in Puis et al. (2005). In order to reduce convergence 
times, the theoretical calculations included only the H and He 
or the H, He and Si model atoms (depending on the strength 
of the Si lines in the spectra and their usefulness as a Toff di- 
agnostic). However in order to incorporate the line blanketing 



from important metals not included in the rate equations an ap- 
proximate treatment was invoked (see Puis et al. for details and 
comparisons with a more detailed technique used by the model 
atmosphere code CMFGEN). The analysis method described be- 
low follows closely the iterative processes outlined in HeiTero, 
Puis & NajaiTO (2002 1 and references therein. 

For the majority of the O-type stars the effective tempera- 
ture was determined by profile fitting of the neutral and singly 
ionized helium spectra, with more emphasis given to the latter 
For the cooler objects in this subset (i.e. 09 - BO), estimates 
were also obtained using the silicon ionization equilibrium (Si 
Ill/iv), with good agreement being obtained between the two 
approaches. Profile fitting to the wings of the hydrogen lines 
was used to determine the surface gravity. Greater weight was 
given to the higher order Balmer lines where the wind has less 
influence on the core of the lines. The wind properties were de- 
scribed by the mass-loss rate, the /3-parameter and the terminal 
velocity, the latter was adopted from the spectral-type - termi- 
nal velocity scale introduced by Kudritzki & Puis (I2000> . The 
/3-parameter was initially assumed to be 0.8 and iterated along 
with the mass-loss rate to obtain a good fit to the H^ profile. 
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Fig. 2. An example of the chi-square minimization procedure 
used to estimate the projected rotational velocity for the star 
4755-030. 



4755_030 




150 200 250 300 350 400 450 
^sini [km/s] 



3.2. Projected rotational velocities 

Stellar rotational velocities for the B-type stars were deter- 
mined by comparing observed line profiles to theoretical coun- 
terparts that had been convolved with a rotational broadening 
function (see, for example, Grav l 19921 . Although the hydrogen 
lines are well observed throughout our sample, they are intrin- 
sically broad, making it difficult to disentangle effects due to 
rotation from those due to the atmospheric parameters, whilst 
there are no metal lines that are usable across all the observed 
B spectral sub-types. The helium spectrum was a suitable com- 
promise between being well observed whilst not being too in- 
trinsically broad. In particular the He I 4026A line was rela- 
tively strong and thus measurable in most of the program stars. 
Its theoretical line shape also appeared to be well calculated 
when its forbidden components were taken into account. In ad- 
dition, the line having been already used for the temperature 
estimation provided an additional consistency in the analysis. 

For each program star, a theoretical spectrum was selected 
at a matching, or the nearest, point in our grid of non-LTE 
model atmospheres to the estimated atmospheric parameters 
with a microturbulence of 5 km s~^ again being adopted. To 
allow for any difference in line strength between the model and 
the observed spectra, the ratio of their equivalent widths (which 
was generally near to unity) was used to scale the model spec- 
trum. We note that a broader continuum range of ±10A from 
the line centre was adopted than that used for the estimation of 
the atmospheric parameters. This was to ensure the inclusion 
of the wings of the profile which were important when esti- 
mating the rotational velocity. The scaled model spectrum was 
then convolved with various projected rotational velocity val- 
ues, and the best fit to the observation, was obtained through a 
chi-square minimization test (see Fig.|2lfor an example). This 
procedure is a simplified version of that used by Ryans et al. 
(2002 ) to estimate the contribution of rotation and macrotur- 
bulence in the broadening of the spectra of early-type super- 
giants. Tests reported by Ryans et al. indicate that the scal- 
ing of the equivalent widths and the values adopted for the at- 



mospheric parameters (including the microturbulence) are not 
likely to lead to significant errors for stars with significant ro- 
tational broadening. Fig. |3] shows examples of the optimal fits 
obtained when estimating the projected rotational velocity for 
two stars with relatively low and high projected rotational ve- 
locities, with all the estimates listed in Tables to |9] having 
been rounded to the nearest 5 km s^^. 

We note that the estimates for the very narrow-lined objects 
may be unreliable due to the intrinsic width of the He I diffuse 
lines. Hence for all stars where we deduced a v sin i of less 
than 50 km s~^, we have undertaken a re-analysis of the stel- 
lar spectra. In these cases the Si III line at 4552Ain the hotter 
stars and the Mgll line at 4481 A in the cooler objects were 
considered. A similar methodology as for the helium lines was 
adopted with the theoretical profiles again being taken from 
our non-LTE grid with a micro turbulent velocity of 5 km s^^ 
being adopted and the doublet splitting for the Mg II line be- 
ing explicitly included. In most cases the values obtained from 
the helium and metal lines were similar but in some cases they 
differed by up to 20 km s^^. For these narrow lined stars the 
values listed in Tables0to|9]are those estimated from the metal 
lines. 

Instrumental broadening has not been included in our anal- 
ysis and we have investigated how this might affect our esti- 
mates as follows. Representative theoretical profiles have been 
convoluted with either just a rotational broadening function or 
with a rotational broadening function and a Gaussian profiles 
to represent instrumental broadening (with a FWHM equiva- 
lent to the inverse of the spectral resolution). For the FLAMES 
data, the profiles only become significantly different for pro- 
jected rotation velocities less than 20 kms~^, with a corre- 
spondingly smaller value for the higher resolution FEROS data. 
Such small estimates of the projected rotational velocity have 
only been found for two targets in our sample and we have not 
attempted to correct these two estimates as they will be sub- 
ject to additional uncertainties due to the choice of the intrinsic 
profiles and in particular the adopted microturbulence. 

For the O-type stars in NGC 6611, the projected rota- 
tional velocities were determined simultaneously with the at- 
mospheric parameters (see Sect. I3.1.3I I using a method similar 
to the B-type stars. For these stars the N III lines at 45 10-45 1 8 A 
and the Mg II at 4481 A line were considered. 

3.3. Correction of Atmospheric Parameters 

Fig- S shows the distribution of our initial estimates of the at- 
mospheric parameters of stars observed in NGC 3293, deduced 
from the hydrogen and helium line strengths as discussed in 
Sec. 13.1.11 The distribution of surface gravities is surprising in 
that there appears to be a significant number of low gravity ob- 
jects. Given the age of this cluster, it would be expected that 
main sequence stars would make up the most of the sample for 
this range of spectral types. In particular forty percent of the 
stars appear to have surface gravities less than 3.70 dex, which 
is closer to that expected for giants or supergiants than to main 
sequence stars. We suspect this may have arisen from the use 
of a restricted wavelength region when measuring the equiva- 
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Table 2. Stars that have been excluded from the analysis. 
Objects later than B8 are not listed and can be identified from 
Paper I. 





3293 


4755 


6611 


Binaries 




24 


7, 13, 14, 30, 68 


Emission 


11 


14, 18 


10, 22, 28 


Others 


72, 121 







lent widths of the hydrogen and helium lines. Stellar rotation 
broadens line profiles, and when a line is normalized at a fixed 
continuum position, the measured EW could be different from 
its intrinsic value, with the discrepancy increasing with the pro- 
jected rotational velocity. As we have compared such observa- 
tional equivalent widths with values measured from theoreti- 
cal spectra with no rotational broadening, the estimated atmo- 
spheric parameters might be incorrect. 

To investigate this effect, we have constructed a grid of 
correction factors for different projected rotational velocities. 
Theoretical spectra have been convolved with rotational broad- 
ening functions and then the equivalent widths of hydrogen and 
helium lines were measured. These have been used to deduce 
atmospheric parameters from our calibration to estimate how 
much they deviate from the actual model parameters. This de- 
viation is indicative of the effects of rotation on our adopted 
methodology, and has been used to correct our initial estimates 
of the atmospheric parameters. For example. Fig. |5] illustrates 
the corrections for a projected rotational velocity of vsini = 
300 km s~^ - theoretical spectra taken from grid points are 
marked with a diamond (o) and after convolution, they im- 
ply the set of atmospheric parameters marked with a plus sign 
{+). The pairs of numbers next to the latter are the derived at- 
mospheric parameters, with the amount of deviation presented 
in parenthesis. Correction factors were also calculated for pro- 
jected rotational velocities from to 360 km s^^, covering all 
the V sin i values of our program stars. The atmospheric param- 
eters of the program stars were then corrected through a series 
of linear interpolations between the nearest grid points in Toff, 
log g and v sin i and these corrected values are summarized in 
Tables0-|9] Tests show that use of these corrected atmospheric 
parameters did not lead to any significant changes in the values 
previously deduced for the projected rotational velocities. 

3.4. Comments on individual stars 

Some stars have been excluded from the analysis as the mea- 
surements of their hydrogen and helium equivalent widths were 
either not possible or unreliable. They largely fall into three 
categories, viz. targets later than B8, binaries or stars with very 
strong emission feature in their Balmer line profiles and are 
summarized in Table |2] This leaves 92 out of 126 stars ob- 
served in NGC 3293, with the numbers for the other clusters 
being 87/108 (NGC 4755) and 44/85 (NGC 6611). Some com- 
ments on individual stars are given below: 



3293-001, -002 and 4755-004: the effective temperature of 
these supergiants were estimated utilizing the silicon ionization 
balance. Their microturbulence has been determined to be 20, 
14 and 15 km^^, respectively, by minimizing the slope in the 
abundance estimates from the Si III lines versus their equivalent 
widths. 

3293-027, -029 and -032: the spectral types of these stars 
remain uncertain due to their very broad line profiles, and they 
were classified as BO. 5- 1.5 in Paper I. We assume that they are 
B 1 .0 and adopt a 'calibrated' effective temperature of 23 000 K 
(TableQ. 

3293-034: The helium spectrum implies that this B2 Illh 
star has an effective temperature of 26 000 K, much higher than 
the limit of validity for this method and incompatible with its 
spectral type. We suspect that these inconsistencies arise from 
its spectral peculiarities. 

6611-010 & -022 These Be stars have no useful effective 
temperature indicators, making it difficult to estimate their at- 
mospheric parameters. Theoretical profiles were fitted to H7 
line to produces loci of possible atmospheric parameters. For 
representative parameter pairs, usini values were estimated 
and their average is presented in Table|9l 

3.5. Stellar masses and luminosities 

In order to determine stellar luminosities, it is important that we 
have reliable cluster distances. Several studies have estimated 
these distances and some of the more recent are summarized 
in Table|3] We have chosen to list results from wide-field CCD 
and spectroscopic surveys, which supersede the older photo- 
electric methods based on less reliable samples. 

We have also independently determined a distance to each 
cluster from a spectroscopic parallax method, applied to the 
main-sequence targets in our clusters. Firstly intrinsic colours 
(B — V)o were calculated using standard Johnson BV filter 
curves and flux calibrated TLUSTY spectra for a range of ef- 
fective temperatures and two gravities. As can be seen from 
Table |5]there are relatively small differences between the val- 
ues for log 5 = 3.0 and 4.0 dex . Hence we adopted the former 
for all stars with log 5 < 3.5dex, and the latter for the higher 
gravity objects. A standard Galactic reddening law of A„ = 
3.1 E(B - V) was adopted for NGC 3293 and NGC 4755 as 
this has been found previously to be a suitable value (Baume et 
ar2003V However there is evidence that the reddening towards 
NGC 661 1 is anomalous and a value of Ay = 3.75 E(B - V) 
has been found to best match the optical and infra-red magni- 
tudes of the cluster members when combined with the colour 
excesses from stars with known spectral types (Hillenbrand 
IT9931 . 

From the non-rotating evolutionary tracks of Schaller et 
al. (.1992 J and Meynet et al. ( 1994.) we have inferred surface 
gravities and then by interpolating in the log g-Tcs plane, es- 
timated stellar luminosities. These have been combined with 
bolometric corrections and the E{B — V) values to determine 
spectroscopic distances for each target. For log T^s >4.45, the 
bolometric corrections of Vacca et al. ( 2004 1 were used, while 
for lower effective temperatures, those of Balona Cl994j were 
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Cluster 


Method 


Distance 


Age 






kpc 


Myr 


NGC 3293 








Frey hammer et al. (.2005 ) 


Eclipsing Binary 


2.8±0.3 


10-13 


Baume et al. <2003l 


UBVRI MS fitting 


2.8±0.3 


6.5-10 


Shobbrook (m3\ 


uvhyP sequence fitting 


2.5±0.1 




This paper (33 stars) 


Spectroscopic parallax 


2.9±0.1 


10-20 


Mean of all results 




2.8±0.2 


10 


Mean and SD of Ay 




0.80±0.29 





NGC 4755 

Sanner et al. (T995t 

Balona & Koen ^1994^ 
Sagar & Cannon ( 1995 I 
This paper (39 stars) 
Mean of all results 
Mean and SD of Ay 



UBVRI MS fitting 
uvbyP MS fitting 
UBVRI MS fitting 
Spectroscopic parallax 



2.1±0.2 8-12 
2.0±0.2 

2.1±0.2 10 

2.3±0.1 10-15 

2.1±0.2 10 
1.15±0.15 



NGC 6611 

Hillenbrand et al. ^T993^ 
Belikovetal. (19991 
This paper (24 stars) 
Mean of all results 
Mean and SD of Av 



Spectroscopic parallax 2.0±0.2 1-3 
(751/ MS fitting 2.1±0.2 

Spectroscopic parallax 1.8±0. 1 2-4 

2.0±0.2 2 
2.81±0.88 



adopted. As discussed above, the determination of the surface 
gravity in our fastest rotating targets is subject to some uncer- 
tainty. Hence we restricted our sample of stars to those, hav- 
ing 4.5 < \ogg < 3.95. Effectively we are estimating cluster 
distances by applying a hybrid spectroscopic and evolutionary 
parallax method for stars which we are confident are lying close 
to the main-sequence and have well determined surface grav- 
ities. The distance estimates are listed in Table|3] and are in 
good agreement with those found previously from other meth- 
ods. We have adopted a mean value of these estimates for use 
in this paper 

Using these distances, we have re-determined the luminos- 
ity of each star from their F-band magnitude, a bolometric cor- 
rection (as discussed above) and their E{B ~ V) value. This 
allowed stars to be placed on a HR-diagram, and by interpola- 
tion between the Geneva evolutionary tracks led to estimates of 
their evolutionary masses, which are summarized in Tables 0- 

m 

3.6. Cluster membership 

In order to investigate whether an individual star is a clus- 
ter member we have used two criteria. One is radial velocity 
(Usted in TablesEHSl and taken directly from Paper I) and the 
other is distance from the cluster centre. In Table |3 we sum- 
marize the mean radial velocity of the clusters (again as re- 



ported in Paper I), together with the adopted cluster centres and 
radii (NGC 3293: Baume et al.lOOT, NGC 4755: Lvnga [T987l 
NGC 66 11: Belikov et al. .1999i) . The use of radial velocities 
alone is not a reliable way to determine cluster membership as 
stars with discrepant (but apparently non-variable) radial ve- 
locities could be binaries with periods significantly longer than 
our sampling frequency (which for NGC 3293 and NGC 4755 
is less than one night). Hence we have only identified a star as 
a non-member if its radial velocity is more than 2cr away from 
that of the mean value and its distance from the adopted clus- 
ter centre (designated r' in Tables is more than twice the 
cluster radius. Such objects have been identified by a f and we 
note that they are only a small fraction of our total sample. In 
turn this may imply that our rejection criteria may not be suf- 
ficiently stringent but we have found that changing our criteria 
(e.g. to radial distances of less than one cluster radius) had little 
effect on our sample. 

4. Discussion 

4. 1. Atmospheric parameters 

We have used a variety of methods based on non-LTE model 
atmosphere techniques to constrain the stellar atmospheric pa- 
rameters. As discussed in Sect. 13.1.31 the use of different ef- 
fective temperature criteria for the O and BO-type stars led 
to excellent agreement. Additionally previous analyzes of the 
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Table 4. The mean cluster radial velocities (v^), cluster centres 
and radii (see text for details) 



Cluster Vr Centre Radius 

km s"^ RA Dec 

NGC3293 -12±5 lO'' 35" 49.3" -58° 13' 28" 4.1' 

NGC 4755 -20±5 lO'' 53" 39.0" -60° 21' 39" 5.0' 

NGC 6611 +10±8 18'M8'" 40.0" -13° 47' 06" 4.8' 



Table 6. The mean projected rotational velocities for our higher 
gravity targets (log g > 3.9; designated 'High gravity') and our 
apparently evolved targets excluding known supergiants (log g 
< 3.9; designated 'Low gravity'). The number of stars in each 
sample is given by n. 



Cluster 


High gravity 


n 


Low gravity 


n 




km s^^ 


km s^^ 






NGC 3293 


184 ± 99 


47 


196 ± 95 


41 


NGC 4755 


173 ± 71 


47 


211 ±96 


38 


NGC 6611 


132 ± 75 


22 


187 ± 73 


6 


All targets 


170 ± 86 


116 


202 ± 94 


85 



considered. This result appears anomalous as if the lower grav- 
ity targets had evolved from the main sequence, their rotational 
velocities might have been expected to decrease due both to the 
increase in their moment of inertia and to the loss of angular 
momentum from mass loss. Rather this trend may arise from 
the estimation of surface gravities for rapidly rotating stars and 
we have identified two possible causes. Firstly, the correction 
procedure for the atmospheric parameters discussed in Sect. 
I3.3l mav not have been successful. To test this hypothesis, we 
have used the atmospheric parameters and projected rotational 
velocities given in Tables0j9]to generate theoretical hydrogen 
line profiles. In effectively all cases these were in good agree- 
ment with our observed line profiles and hence we believe that 
our correction procedure is not the main cause of this discrep- 
ancy. 

A second explanation is that for stars with large rotational 
velocities, the effective surface gravity will vary over the stel- 
lar surface due to the centrifugal force increasing as one moves 
from the pole to the equator. To a first order of approxima- 
tion, our methodology will measure the mean effective sur- 
face gravity, which will be smaller than that of an analagous 
slowly rotating object. Hence the anti-correlation between pro- 
jected rotational velocity and estimated surface gravity could 
be a manifestation of such an effect. For example, Howarth & 
Smith (2001 1 found a variation of 0.4 dex between the effec- 
tive surface gravity at the pole and the equator for three rapidly 
rotating 09-09.5 type stars. 

To test this hypothesis we have computed models and syn- 
thetic spectra for a star with an B spectral type (corresponding 
to rcfr=20000; log (7=4.0 at the stellar poles) using the method- 
ology discussed by Townsend (1997i which includes the ef- 
fects of rotation on the pole to equator temperature structure 
and on the stellar radius. The synthetic input spectra were based 
on LTE atmospheres, calculated with ATLAS9 (Kurucz|2005} 
and LTE line formation, which should be adequate for assess- 
ing the effects of stellar rotation. Further details can be found 
in Kauferet al. (2006). 

For models with different rotational velocities, we have 
measured the equivalent widths of the features used to estimate 
atmospheric parameters such as the He I line at 4026A and the 
hydrogen B aimer lines. The results are difficult to interpret but 
the effects of rotation on these spectral features appear negli- 
gible for equatorial rotational velocities less than 250 km s~^. 



spectra of B-type supergiants using the TLUSTY grid adopted 
here (Dufton et al. I2005> and the code fastwind (Trundle et 
al. '20051 have also yielded encouraging agreement. Hence we 
believe that in general our estimates of effective temperature 
should have an accuracy of typically 5% and should be gener- 
ally more reliable than those based on the calibrations of broad 
band colours. Exceptions may occur for the small number of 
stars that are peculiar (see Sect. 13. 4> or where we have adopted 
effective temperatures based on spectral types. 

For our gravity estimates, the situation was complicated by 
the corrections required for targets with large projected rota- 
tional velocities (see Sect. 13. 3> . Whilst the corrections for the 
effective temperature estimate are relative small (typically less 
than 1 000 K), those for the surface gravity can be significant 
(see Fig. |5}. To investigate the reliability of these corrections 
we have considered how the gravity estimates vary with rota- 
tional broadening. Our procedure has been to divide our targets 
into two bins — those that appear to be near to the main se- 
quence (log g > 3.9 dex) and those that may have evolved. For 
the latter, we have excluded five supergiants (classified as la or 
lb in Paper I), as these targets are clearly evolved. Additionally 
we have excluded the O-type stars for NGC 66 11 as these were 
analyzed separately from the B-type stellar sample and em- 
ployed detailed fitting of the hydrogen line profiles. In Table 
|6l we summarize the mean and standard deviations for the pro- 
jected rotational velocities of our two samples for each individ- 
ual cluster and for all our targets. 

For all three clusters (although the sample sizes for 
NGC 66 11 are small), there is a systematic trend in the sense 
that the low gravity objects have the higher mean projected ro- 
tational velocity and this is confirmed when all the targets are 

Table 5. Johnson UBV colours deduced from the tlusty 
spectra 



log g = 4.0 log g = 3.0 



Tcff 


B~V 


U-V 


B~V 


u ~v 


13500 


-0.102 


-0.402 


-0.119 


-0.478 


15000 


-0.125 


-0.500 


-0.140 


-0.571 


18000 


-0.161 


-0.649 


-0.174 


-0.724 


20000 


-0.181 


-0.732 


-0.188 


-0.813 


23000 


-0.203 


-0.840 


-0.198 


-0.917 


25000 


-0.212 


-0.899 


-0.209 


-0.958 


27500 


-0.228 


-0.947 


-0.209 


-0.989 


30000 


-0.244 


-0.984 






35000 


-0.258 


-1.035 
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For faster rotators, the equivalent width of the He I remains ef- 
fectively unchanged but the strength of the hydrogen lines de- 
creases. However the changes are relatively small with those 
for the HS equivalent width being approximately 10-15% for 
equatorial rotational velocities upto 400 km s^^, which corre- 
sponds to 80% of the critical equatorial velocity. In turn this 
would lead to a change in the estimated surface gravity of 0.2 
to 0.3 dex. 

In Sect. 14.51 we estimate the mean equatorial rotational 
velocity for the cluster NGC 3293 and NGC 4755 as approx- 
imately 225-250 km s^^. Assuming a Gaussian distribution of 
velocities, the -i width would be approximately 110 kms^^. 
Hence the number of targets with equatorial velocities greater 
than 250 km s^^ would be 40-50%, whilst those with a rota- 
tional velocity of 400 km s^^or more would be of the order 
of 5-10%. From the above, although high rotational velocities 
may contribute to the spread of gravity estimates found in our 
two older clusters, it is unclear whether it can fully explain 
these variations. 

In summary, we believe that our effective temperature es- 
timates are secure and find no evidence that our methodology 
has introduced systematic errors in our gravity determinations. 
Rather the significant number of relatively low gravity B-type 
objects may reflect in part a decrease in the effective surface 
gravity due to centrifugal acceleration. 

4.2. Hertzsprung-Russell diagrams 

For all three cluster, we have used the luminosities and 
effective temperatures listed in Tables EHU ^ construct 
Hertzsprung-Russell (HR) diagrams and these are presented in 
Fig.|6l We have used different symbols to distinguish between 
those targets with estimated logarithmic surface gravities less 
than and greater than 3.9 dex. These are effectively the same as 
the samples discussed in Sect. l4.1l but with the supergiants now 
being included in the former. Also shown are the evolutionary 
tracks and isochrones of Schaller et al. ( 1992 1 and Meynet et 
al. J1994 >. taken from the database of the Geneve group (see 
[http://obswww.unige.ch/~mowlavi/evol/stevjdatabase.html 
also available from the Centre de Donnes Astronomiques 
de Strasbourg, http://cdsweb.u -strasbg.fr/CDS.htmli ) for a 
metallicity, Z = 0.02. 

In general there is a good qualitative agreement between 
the observed and theoretical HR diagrams although the ob- 
served main sequences appear to be systematically more lumi- 
nous than predicted. For the late-B type stars, where the choice 
of isochrone is not important, this discrepancy would appear 
to be approximately, Alog L~ 0.4, for the older clusters. The 
discrepancy may be larger for NGC 66 11 but in this case the 
sample size is relatively small. This systematic difference in 
luminosity would correspond to an error in our adopted dis- 
tances of approximately 60%, which although possible is sig- 
nificantly larger than the range of values summarized in Table 
|3l Slesnick et al. (1992) have undertaken an extensive pho- 
tometric and spectroscopic study of the double cluster h and 
X Persei, which is similar in age to those of NGC 3293 and 
NGC 4755. Inspection of the HR diagrams (see their Fig. 6 



and 7) imply good agreement between the observed late-B type 
stars and the theoretical predictions of Schaller et al. J1992> . 
However, recently, Strom et al ( 12005 1 have undertaken a study 
of projected stellar rotational velocities in h and x Persei and 
used the results of Slesnick et al. For this smaller sample, there 
is again evidence of a discrepancy (see Fig. 1 of Strom et al.) 
similar to that found here for the later B-type spectral types. 

In our HR diagrams, we have distinguished between appar- 
ently low and high gravity objects as discussed above. For our 
two older cluster, the effects of stars evolving from the main 
sequence can be seen for our early-B type objects. However 
for the later B-spectral types, there is no clear segregation of 
the targets as would be expected if the lower gravity objects 
had evolved from the main sequence. We have further investi- 
gated this by refining the analysis discussed in Sect. l4.1l to only 
consider targets with effective temperatures less than 20 000 K, 
where evolutionary effects should be negligible. For all our tar- 
gets, the mean projected rotational velocities are 178±72 and 
217±92km s^^for the high and low gravity stars respectively. 
The former is similar to that listed in Table |6l for the corre- 
sponding sample, whilst the latter is larger due to the exclusion 
of hotter evolved objects. The position of the cool lower gravity 
objects in the HR diagrams and their high mean rotational ve- 
locities both support the hypothesis that our gravity estimates 
are reliable (in the sense that they are measuring the mean sur- 
face gravity) and that their range may reflect the different con- 
tributions of the centrifugal acceleration. An important conse- 
quence of this is that the lower effective gravity has led to many 
of these objects being classified as giants (see Tables0]9} and 
hence calibrations based on spectral type of, for example, the 
luminosity would led to systematic errors. 

From the HR diagrams it was also possible to estimate clus- 
ter ages (assuming coeval star formation). For all three clusters 
the fitting of the isochrones was problematic as discussed be- 
low. 

NGC 6611: As discussed above, the observed main se- 
quence (which for this young cluster is still discernible at BO 
spectral type) appears systematically brighter than predicted. 
There is evidence that the extinction towards this cluster is 
anomalous (see, for example, Chini & Wargau, 1990, and ref- 
erences therein) and indeed we adopted a different value for 
the ratio of the total extinction to the reddening, R (see Sect. 
13. 5> than for the other clusters. Hence some of the discrepancy 
may arise from uncertainties in the interstellar extinction. We 
note that the HR diagram of Hillenbrand et al. (1993' Fig. 8) 
implies a similar discrepancy. These authors also considered 
isochrones for pre-main sequence evolutionary tracks but for 
our targets with effective temperature greater than 20 000 K the 
timescales for reaching the main sequence would appear too 
small to explain this discrepancy - for example 300 000 years 
at an effective temperature of 20 OOOK. However Hillenbrand et 
al. also found evidence that some intermediate mass stars (3-8 
solar masses) could be as young as 250000 years. If such ages 
were also present amongst the more massive stars, this would 
help explain the enhanced luminosity of our objects. We find a 
best fit isochrone for this cluster of approximately 2 an , with 
some evidence for age spread of 1 Myr, consistent with the re- 
sults of Hillenbrand et al. 
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Fig. 6. Hertzspmng-Russell diagrams for the three clusters, together with evolutionary tracks and isochrones taken from the 
Geneva database. Stars with estimated gravities of less than 3.9 dex are shown as open circles with higher gravity stars represented 
by closed circles. 
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NGC 3293 and NGC 4755: Given the similarity in their 
HR diagrams, these clusters are best considered together. As 
discussed above, their observed main sequences (for effective 
temperatures of less than 20 000 K) appear to be more lumi- 
nous than predicted by evolutionary models. In these cases, 
this would not appear to be due to anomalous extinction. 
Also, given the larger cluster ages, it would appear unlikely 
that we are observing stars evolving to the main sequence, 
although Sagar & Cannon (il995t identified lower mass pre- 
main sequence objects with ages between 3-10 million years 
in NGC 4755. Inspection of their Fig. 7 implies that Sagar & 
Cannon also found the B-type main sequence to be overlu- 
minous compared with the isochrones of Maeder & Maynet 
<1991> . By contrast. Fig. 7 of Baume et al. ( 2003 ) implies good 
agreement between the observed HR diagram of NGC 3293 
and the theoretical models of Girardi et al. (2000 1. However 
this comparison is complicated by Baume et al. extending the 
envelope of the predicted B-type main sequence upwards by 
0.75 magnitudes corresponding to binary systems with two 
similar stars. Additionally the published isochrones of Girardi 
et al. have a minimum age larger than those quoted by Baume 
et al., who presumably had access to additional unpublished 
information further complicating any comparisons. A binary 
population could also help explain our discrepancy but given 
our use of high dispersion spectroscopy to estimate effective 
temperatures (compared with Baume et al. who used (B-V)o), 
we would expect any such contamination to be small. From 
the position of the main sequence turnoff, both clusters ap- 
pear to have ages between 10 and 20 milUon years. The the- 
oretical isochrones fit the giant population reasonably well in 
NGC 4755 but in NGC 3293, there are a group of stars with 
an implied evolutionary age greater than 20 million years. 
Additionally both clusters contain supergiants (the highest lu- 
minosity stars in Fig. |6j, which if they have evolved as single 
stars appear to have evolutionary lifetimes of less than 10 mil- 
lion years. Indeed the lower age estimate obtained by Baume et 
al. may be at least partly influenced by these targets. 



In summary, for all three clusters, the agreement between 
the observed HR diagrams and the theoretical predictions is 
rather disappointing. This is particularly the case as our effec- 
tive temperature estimates should be more reliable than other 
studies that have used photometric colours. Indeed for all three 
clusters it would appear that there is either non-coeval star for- 
mation or the evolution of individual targets is being influenced 
by other factors. The former is a plausible explanation for 
NGC 66 1 1 given its youth and the detection by Hillenbrand et 
al. (1993) of very young pre-main sequence intermediate mass 
objects. Although the bulk of the star formation for NGC 3293 
and NGC 4755 is consistent with ages of 15Myrs, the most 
massive stars in the cluters have smaller evolutionary ages of 
approximately 5 Myrs. This would imply a very extended pe- 
riod of star formation and hence these discrepancies may be 
due, at least in part, to other factors such as binarity. 



4.3. Initial Mass Functions 

The atmospheric parameters and masses determined in the pre- 
vious sections allow the mass functions of the three clusters to 
be re-evaluated with the most extensive data and information 
to date. As discussed in Paper 1, we are not spectroscopically 
"complete" in any of the clusters, in the sense that we were not 
able to place fibres on all stars within a given colour-magnitude 
range. However in the case of NGC 3293 and NGC 4755 the 
reddening towards the cluster appears fairly homogeneous and 
hence we can estimate the effective temperatures and luminosi- 
ties from the stellar BV photometry. In the case of NGC 661 1, 
which shows a scattered colour-magnitude diagram, such a pro- 
cedure would not be reliable and we have chosen to increment 
the Hillenbrand et al. fl993'( study with our additional spec- 
troscopy in order to improve the statistical significance of the 
results. 

4.3.1. NGC 3293 

We have selected all the stars from the colour magnitude di- 
agram (see Paper I, Fig. 3) for which we do not have spec- 
troscopic data and which were in the range V < 13.3 and 
{B — V)o < (where {B — V)o is the intrinsic colour esti- 
mated using the data in Table |3}- Effective temperatures were 
estimated by fitting a cubic to the Tcs~{B — V)q results in 
Table |5] (assuming that log 5 = 4.0). We then calculated bol- 
metric luminosities and masses using the method employed in 
Section lTsl Thirty six stars were selected, all with masses in 
the range 3-10Mq. There is also one confirmed red supergiant 
cluster member of NGC 3293 which must be included. Feast 
( tl958) gives an MO lab classification with A/y = — 5.9 and us- 
ing the bolometric corrections of Elias et al. J1985> . we have es- 
timated a mass from the Geneva evolutionary tracks of 15M0. 

To determine the IMF, we follow the method of Hillenbrand 
et al. (1993 1 and Massey et al. ([19951, so that our results can 
be directly compared to those studies. Stars were divided into 
mass bins, whose size was chosen to give comparable numbers 
of objects as far as possible. The number of stars in each bin 
was normalised to unit logarithmic mass interval and to unit 
area (kpc^) and the Scalo (1986) notation was adopted. The 
area sampled was taken from the cluster radius (we used twice 
the quoted radius, to be consistent with the membership criteria 
applied above). 

In Fig.0we plot the present day mass function (PDMF), 
which given the relatively small age of NGC 3293 and the fact 
that the majority of stars are on the main-sequence, is effec- 
tively the IMF. A linear least squares fit was performed, with 
each point being weighted by vW, where N is the number of 
points in the bin. We find a slope ofF = — 1.5±0.2 over the 
range 2.8-4OM0. Baume et al. ( l2003t estimate F = -1.2 for 
the mass range 1.4-45M0, but find a steeper slope of F = —1.6 
for the higher mass stars in the range 8-45M0. Our results 
imply that this steeper value for the high mass objects ex- 
tends down to ^3M0 and are consistent with the range of 
IMFs found by Massey et al. J 19951 1 for young clusters in the 
Northern Milky Way. 
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Fig. 9. Cumulative probability function for all targets in our two older clusters and in h and x Persei and for those in the different 
classes as specified by Strom et al. (2005). Group 1 - relatively unevolved stars with an effective temperature, Teff < 17000; 
Group 2 - stars with 17000 < Tcfr< 21000; Group 3 - stars that have evolved significantly from the main sequence, Tes> 21000. 
The solid lines are for field stars, the dotted lines are for the h and x Persei sample of Strom et al. J2005> and the dot-dash lines 
are our combined sample for the clusters NGC 3293 and NGC 4755. 
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4.3.2. NGC4755 

The PDMF for NGC 4755, and hence by implication the IMF, 
was estimated in a similar way as that for NGC 3293 . We 
used the same magnitude and colour criteria {V < 13.3 and 
{B — V)o < 0) and identified 51 stars, with masses in the 
range 2-4M0. Again there is one red supergiant that is a cluster 
member, with an M21b spectral type (Feast 1963 ), for which we 
have estimated a mass of 15M0. In Fig.0 we plot the IMF for 
NGC 4755 and estimate a slope of r = -1.8±0. 2, which again 
is consistent with the range reported in Massey et al. (.1995 J and 
similar to that for NGC 3293. 

4.3.3. NGC 6611 

For NGC 66 1 1 , because of the highly variable reddening across 
the field, it was not possible to use a similar method to that for 
NGC 3293 and NGC 4755. From Hillenbrand et al. (>1993j . we 



can identify 22 stars with spectral types earlier than B3, which 
are likely cluster members but for which we do not have spec- 
troscopy. All but one of these are classified as B0.5V-B3V and 
we have adopted effective temperature estimates based on the 
data in Tabled The other target is BD— 13°4912, which was 
classified as B2.5I, and for this we have adopted Toff — 16 500 
(Crowther et al. 2006 1. Masses were then calculated following 
the same procedure as discussed in Section l331 

We have also estimated effective temperatures, luminosities 
and masses for likely cluster members in the Hillenbrand et 
al. sample without spectral types by calculating Q, and using 
the Q — Teff relation of Massey et al. (19_89 l. None of these 
targets had masses above 7Mq and hence it is reasonable to 
assume that the combination of our spectroscopy and that of 
Hillenbrand et al. is virtually complete above 7Mq. In Fig.Q 
we plot the IMF for NGC 66 1 1 and determine a slope of F = 

— 1.5 ± 0.1, which is slightly steeper than the value of F = 

— 1 . 1±0.3 of Hillenbrand et al., but is consistent with a Salpeter 
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type value and again similar to our other clusters and those in 
the Massey et al. (.1995) survey. 

4.4. Projected Rotational Velocities 

In Figure |8l the distribution of projected rotational velocities 
are plotted as histograms for all three clusters. Given the simi- 
larity in their ages and HR diagrams, it is not surprising that 
the distributions for NGC 3293 and NGC 4755 appear simi- 
lar. There appears to be an excess of stars in NGC 4755 with 
projected rotational velocities between 100-150 km s^^but as- 
suming Poisson statistics, this is not statistically significant and 
indeed the mean values for the two clusters are almost iden- 
tical (NGC 3293: 188 kms'^ NGC 4755: 185 kms"^). We 
have calculated cumulative probability functions (CPFs) for 
our combined samples in NGC 3293 and NGC 4755 in order 
to compare our results with the study of Strom et al. (120051 ) 
for the double cluster, h and x Persei, which appears to have a 
similar age (of approximately 13 million years for their nuclear 
regions - Slesnick et al.,'1992'). We follow the methodology of 
Strom et al. by dividing our sample into three groups (see Fig. 
1 of Strom et al.), viz. (1) relatively unevolved stars with an ef- 
fective temperature, Tcff < 17000; (2) stars with 17000 < Tcff< 
21000; (3) stars that have evolved significantly from the main 
sequence, Tes> 21000 K. In Fig.|9l we plot the CPFs for all 
our targets and for each group, together with the CPFs for stars 
in h and x Persei and Strom et al's sample of field stars. Strom 
et al. make two major conclusions about their CPFs. The first 
is that the clusters stars on the whole tend to rotate faster than 
the field stars counterparts. Secondly they find that by dividing 
their sample into these different mass bins they can see a trend 
in which the lower mass stars (4 — OM©) have vsini values 
clearly higher than the field stars but that the higher mass stars 
(9 — I5M0) show a very similar v sin i distribution to their field 
star counterparts. 

Inspection of Fig. |9l shows that the projected rotational 
velocities are systematically higher in the NGC 3293 and 
NGC 4755 clusters than in the field populations. This differ- 
ence has also been found previously in other studies (see, for 
example, Bernacca & Perinotto [T974l Wolff et al. [T982l Keller 
[2004) and has been interpreted in several ways including differ- 
ences in the ages of the field and cluster targets. However our 
results are significantly different to the detailed work by Strom 
et al. Their CPFs for the h and x Persei double cluster and the 
field show a clear differentiation for group (1) targets, decreas- 
ing for groups (2) and (3). Strom et al. discussed this behaviour 
at length and postulated that it reflected the star formation envi- 
ronment rather than evolutionary effects during the stellar life- 
times. Their conclusion that evolutionary effects were unlikely 
was based on the fact that the mean v sin i of the high mass 
stars (group 3) was nearly a factor of two higher than that for 
the lowest mass objects (group 1). Meynet & Maeder (2000) 
and Heger & Langer (2000) calculations of the rotational ve- 
locity evolution of 12 — ISM© stars on the main-sequence 
suggests a decrease in velocity of between 25 — 30% within 
10 — 15yrs. Hence Strom et al. speculated that the differences 
reflected the effects of higher accretion rates that may be char- 



acteristic of star formation in these dense clusters. In particular 
they suggested that enhanced aggregation will have three ef- 
fects, viz. (1) higher initial rotation speeds, (2) higher initial 
radii along the stellar birth line (when deuterium shell burn- 
ing commences), resulting in greater spin up as the star evolves 
to the main sequence and (3) that the higher initial birthline 
radii are especially pronounced for mid- to late-B-type stars. 
This latter leads naturally to the larger differences with respect 
to field stars observed for their group (1) objects in h and x 
Persei. 

In our combined sample for NGC 3293 and NGC 4755, the 
rotational velocity distribution of our group (1) objects is es- 
sentially the same as that found by Strom et al. However our 
group (2) and (3) objects show similar enhanced projected ro- 
tational velocities and in particular do not tend towards those of 
the field stars. Hence contrary to the results for h and x Persei, 
any effect would appear to be operating throughout our stellar 
sample. Huang & Gies (2005) have surveyed 496 OB stars in 
clusters within the approximate age range 6—73 Myrs, and also 
find fewer slower rotators in the cluster stars than in the field. 
They also present some evidence for the higher mass stars (late- 
O to early-B) spinning slower than the lower mass late B-types, 
and suggest that this can be explained by the spin-down of mas- 
sive stars through angular momentum loss via stellar winds e.g. 
as modelled by Meynet & Maeder ( 2000j and Heger & Langer 
(2000). 

One possible explanation for the different projected rota- 
tional velocity distributions for groups (2) and (3) in our com- 
bined NGC 3293 and NGC 4755 sample and that of Strom et 
al. (2005) for h and x Persei is that the clusters have different 
ages. Hence we have used the data presented by Strom et al. to 
construct an HR diagram similar to those shown in Fig.|6l Our 
procedure was to use their absolute visual magnitudes and the 
bolometric corrections of Balona ( 1994) to deduce stellar lu- 
minosities and the HR diagram for h and x Persei is shown in 
Fig- El together with isochrones again taken from the Geneve 
database. 

The observational HR diagram for h and x Persei is qualita- 
tively similar to those for NGC 3293 and NGC 4755. However 
for the former the agreement between observed and theoretical 
main sequences is better for the later B-type stars. As for our 
clusters, the evolved star populations in h and x Persei imply 
a range of ages with an estimate of approximately 15 million 
years being reasonable, consistent with the age of 12.8±1.0 
million years deduced by Slesnick et al. ( I1992> from pho- 
tometry and moderate resolution spectroscopy. These estimates 
are similar to our adopted ages for NGC 3293 and NGC 4755. 
Strom et al. postulate that their projected rotational velocity 
distributions were present when the stars were formed. Hence 
given that the ages of NGC 3293 and NGC 4755 are compara- 
ble with that of the h and x Persei sample, the differences in 
projected rotational velocity distributions cannot be explained 
by variations in the stellar ages of the samples. 

The mean of the v sin i estimates for all the observed stars 
in our younger cluster NGC 66 11 is 144 km s^^, which is lower 
than that for NGC 3293 and NGC 4755 . In NGC 66 1 1 there are 
wsini measurements of 25 stars with masses less than I5M0 
and 13 stars with masses between I5-6OM0, with mean val- 
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Fig. 10. Hertzsprung-Russell diagrams for the h and x Persei cluster, together with evolutionary tracks and isochrones taken from 
the Geneva database. 



ues of 192 and 1 16km s^ ^respectively. That for the lower mass 
group is in good agreement with those for both NGC 3293 
and NGC 4755, where we should be sampling the same stel- 
lar mass range albeit at a different age. For the massive stars, 
evolutionary models (Meynet & Maeder l2000l Heger & Langer 
12000 1 predict that initial rotational velocities of 300km s^^are 
reduced to between 60-200km s^ ^during the main-sequence 
phase due to angular momentum loss via their winds, with the 
most massive luminous objects showing the most pronounced 
effect. The fact that we find a lower mean value of v sin i for 
the more massive stars supports the existence of this effect. 

We can test this hypothesis in a simple and illustrative way. 
For the 13 stars in the 15-6OM0 range we assume that they had 
an initial rotational velocity of 250km s^^(see Sect. 14. 5t and 
estimate the rotational velocity for each star after 2 Myrs (be- 
ing our estimated age for NGC 661 1) using the mass dependent 
evolution of surface equatorial velocities predicted by Meynet 
& Maeder (fM50^ and Heger & Langer (2000). This gives a 
mean equatorial velocity of 166 kms^^, and hence a mean 
V sin i value of 125km s^^ (assuming that the spin axes are ran- 
domly distributed). This is similar to the observed mean v sin i 
of 116kms^^, which suggests that an initial typical equato- 
rial velocity of about 250km s^^and spin down according to 
the evolutionary models is consistent with that observed. We 
caution that this is an illustrative and simple comparison as it 
assumes a single value for the initial velocity and ignores any 
macroturbulent broadening in the atmospheres of the luminous 
objects. 

The mean v sin i of the stars with masses less than 
15M0in each of the three clusters is consistently around 
190km s^^, which would suggest a mean equatorial velocity of 
240km s^^(also see Section l43i . The initial velocity of these 
objects is then likely to be in the range 250-300kms"\ de- 
pending on the evolution of the rotational velocity over the 



course of 2-15Myrs. At present this is uncertain for models less 
than 12M0, as there are no detailed published models for the 
evolution of the rotational velocity. For the 12-15M0 models, 
the rotational velocity does decrease during the first 1-2 Myrs 
but then remains effectively constant for the remainder of the 
main sequence lifetime. Hence it may be reasonable to compare 
the lower mass NGC 661 1 targets with those in NGC 3293 and 
NGC 4755. 

Recently Wolff et al. f2006^ have investigated the pro- 
jected rotational velocities in a sample of very young stars with 
masses ranging from 0.2 to 50 M©. Their principal aim was to 
investigate the formation mechanism for high mass stars and 
in particular if it differed from that for lower mass stars. For 
the full range of spectral types they found that the ratio of the 
projected rotational velocity to the equatorial breakup velocity 
showed no significant variations. Assuming a random orienta- 
tion of rotational axes, this led to a median value for the ratio 
of the rotational velocity to the equatorial breakup velocity, R, 
of 0.14, with values forB-type stars (8 < M < 25Mq) of 0.13 
and forO-type stars (25 < M < 90Mq) of 0.20. 

A similar analysis can be undertaken for our samples with 
the small age and presence of O-type stars making that for 
NGC 661 1 the most appropriate. We have again assumed a ran- 
dom orientation of rotational axes and restricting our sample to 
stars that appear to be main sequence with log g> 3.9 and ob- 
tain a ratio, R, of 0.22 for 27 targets. For our 7 O-type stars, 
the ratio is 0.16, whilst for the 20 B-type stars the ratio is 0.28. 
Hence our ratio for O-type stars is in good agreement with that 
of Wolff et al. but our B-type value is higher although given 
our relatively small sample size this may not be significant. For 
our combined B-type sample for older clusters NGC 3293 and 
NGC 4755, we obtain a far higher median estimate of 0.61 for 
83 targets. The reason for this higher ratio is unclear, as our 
targets would not be expected to have evolved to higher ro- 
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tational velocities during the clusters' lifetimes. However it is 
consistent with the different mean velocities found for the three 
clusters and discussed above. 

In summary for our older clusters, we confirm that the clus- 
ter stars show higher projected rotational velocities than their 
field star counterparts. However we do not reproduce the re- 
sults of Strom et al. for h and x Persei who found that the 
higher mass stars tend to rotate significantly more slowly than 
the lower mass objects. If their explanation of enhanced ac- 
cretion causes the different behaviour between the cluster and 
field samples, then at least for NGC 3293 and NGC 4755 this 
mechanism is present for all B-type objects. The mean equato- 
rial velocity of the most massive stars (I5-6OM0) in NGC 661 1 
is significantly less than the lower mass stars in this cluster 
and those in NGC 4755 and NGC 3293. This could be due to 
the more massive stars losing significant angular momentum 
through stellar winds during the first 2 Myrs of their lifetime 
but is also consistent with the rotational velocities found by 
Wolff et al. C2006 l for very young O-type stars. 

4.5. The intrinsic rotational velocity distribution 

The sample of stars for which we have reliable v sin i values 
is large enough that we may be able to model the intrinsic ro- 
tational velocity distribution by assuming that the inclination 
angle i is randomly distributed in space. Following the formu- 
lation of Chandrasekhar & Miinch ( 1950 1, which is also dis- 
cussed in Brown & Verschueren ( 1997), we employ the inte- 
gral equation which relates the distribution, f{x), of equatorial 
velocities, x, for a sample of stars to the corresponding distri- 
bution of projected rotational velocities (f){v sin i). 

(j){v sini) = V sini / rl.T (1) 

Jvsini x^x^ - (vsmiy 

There are two J{x) distributions for which this equation 
can be integrated analytically. The first is a 5-function and the 
other is a uniform distribution in which all values of x occur 
with equal probability upto a maximum velocity, v^ax- These 
analytic solutions are given in Chandrasekhar & Miinch ( 1950 ). 
We have also considered as an illustrative example a Gaussian 
function for /(.x) and integrated the equation numerically in 
incremental steps of vsmi. Mokiem et al. ( 2 006i have dis- 
cussed the distribution of rotational velocities for a sample of 
O-type stars in the SMC obtained during our FLAMES survey. 
They show that different choices for the distribution of rota- 
tional velocities, /(x), can lead to similar distributions for the 
projected rotational velocities, ^(usini) (see Fig. 10 of their 
paper). Hence we stress that our adoption of a Gaussian distri- 
bution should only be considered illustrative and does not rule 
out the possibility that the actual distribution may be different. 

In order to increase the statistical significance of the model 
fitting we have combined the projected rotational velocity dis- 
ti-ibutions for the samples in NGC 4755 and NGC 3293. These 
two clusters are very similar in age, mass function, metallic- 
ity and total mass as discussed above. Hence we should have 
a sample of stars born in similar conditions and of a compara- 
ble age. This leads to a combined sample of 178 stars and the 



histogram of observed (/)(t;sini) along with the correspond- 
ing CPFs can be compared with our different predictions for 
the projected rotational velocity distribution. We can confi- 
dently rule out the two analytical solutions of the (5-function, 
or the continuous distribution as they provide clearly inconsis- 
tent fits with the observed histogram and CPF as illustrated in 
Fig [n In Fig-El we show the predictions for an intrinsic ro- 
tational velocity distribution, /(x), which is a Gaussian func- 
tion with a peak of 225 km s^^and a \ half width, cr, of 110 
km s^^(corresponding to a full-width-half-maximum, FWHM, 
of approximately 180 km s~^). This is the best fit obtainable 
using a Gaussian profile, both from a subjective comparison 
and a test for a grid of possible solutions. A poorer subjec- 
tive fit (with significantly higher x^) is found when the peak of 
the Gaussian profile is changed by ±25km s^^or the width, a, 
by ±20km s^^. As noted in Section B31 there are 7 stars which 
have evolutionary masses of between 7-4OM0, and which have 
evolved significantly from the main-sequence. When these are 
removed from the sample, we get a slightly higher value of 
250km s^^for the peak of the distribution and a similar width. 
Hence we believe that this is a reasonable solution for the rota- 
tional velocity distribution of stars in the range 3-15Mo at an 
age of approximately 15 Myrs in these two clusters. However 
with our methodology we cannot rule other distributions of a 
similar nature. For example Mokiem et al. (2006) show that a 
Gaussian function would give a similar model distribution as a 
rectangular function of similar width and that a very large num- 
ber of observed stars would be required to definitively distin- 
guish between the two input distributions. Although we cannot 
distinguish between such model distributions, or indeed other 
distributions of similar form, we can conclude that the under- 
lying distribution must be a broad function with a centre and 
width similar to our best fit Gaussian. Additionally we were not 
able to exclude possible binaries from this sample. For example 
Brown & Verschueren ( 1997 ) and Huang & Gies (2005 ) found 
that early-type stars with radial velocity variations on the whole 
rotated more slowly than the likely single star sample. Our ob- 
servational time sampling is not extensive enough to determine 
likely binary fractions, although we should be able to identify 
binaries in our LMC and SMC clusters (Evans et al. 2006). 

In Fig. El ^ similar analysis is shown for NGC 6611. Only 
forty four stars were included, which reduces the significance 
of the fit considerably. A model with a peak of 175km s^^and 
a-=100km s~^(FWHM~ 165 km s"^) appears to be the best fit 
for a Gaussian distribution of rotational velocities. The mean 
rotational velocity is somewhat lower than that for the com- 
bined data for the older clusters, and we note that previous 
studies on the rotational velocities of O-type stars in the Galaxy 
(Howarth et al. 119971 Penny I1996> have indicated that ro- 
tational velocity distributions are peaked around 100km s^^. 
Indeed an analogous Gaussian fit to the unevolved sample of 
Penny (1996) gives a peak of 155km s^^and (7=100kms^^. 
Hence there is a suggestion that the more massive stars (15- 
6OM0) in the Galaxy have typically smaller rotational veloc- 
ity distributions than their lower mass counterparts (3-15M0). 
As discussed in Section l4!4l using a simple illustrative example, 
this can quahtatively be explained by the models of massive 
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rotating stars losing significant angular momentum through 
strong stellar winds (Meynet & Maeder l2003l Heger 8c Langer 

4.6. Comparison with stellar evolutionary models 
including rotation 

The rotational velocity distributions found in the previous sec- 
tion suggest that we should compare our observational HR di- 
agrams to stellar evolutionary calculations which include rota- 
tion. In Fig[0]we reproduce Fig|6l but with evolutionary tracks 
taken from Meynet & Maeder (2003 1 which have an initial ro- 
tational velocity on the ZAMS of 300km s^^(we are grateful 
to Georges Meynet for providing the lower mass tracks 3-7M0 ; 
these are unpublished at present and are computed without core 
overshooting). Meynet & Maeder point out that these models 
would correspond to average equatorial rotational velocities of 
between 180 and 240km s^^during the main sequence phase, 
which is well matched to our preferred solutions for the intrin- 
sic rotational velocity distributions. For the higher mass stars 
we have truncated the tracks before they move bluewards to- 
wards the WR phase in order to aid clarity in the figures. The 
position of the ZAMS in NGC 3293 and NGC 4755 was signifi- 
cantly below the stellar positions for the non-rotating models as 
discussed in Section l4!2l The rotating models do provide some- 
what better agreement, as the ZAMS moves towards lower ef- 
fective temperature. The width of the main-sequence phase also 
increases for a rotating star which may help to explain the scat- 
tering of points above the ZAMS. However one would still ex- 
pect the ZAMS to fit through the middle of the observed distri- 
bution rather than the lower envelope as observed and the cause 
of this is still unclear. If we were to use the rotating models to 
estimate the mass of the stars, we would obtain different val- 
ues to those for the non-rotating models, but not in any simple 
scalable fashion. However the differences are at most 10%, and 
are not likely to influence the PDMF determinations. 

5. Conclusions 

We have analysed high resolution spectroscopy of targets in 
three Galactic clusters to estimate atmospheric parameters and 
projected rotational velocities. These have been used to de- 
rive stellar luminosities and evolutionary masses. Our principle 
conclusions are as follows: 

1 . Atmospheric parameters have been deduced from non-LTE 
model atmosphere calculations and should be more reli- 
able than those inferred from photometry or spectral types. 
Hence the observed HR-diagrams should be most the reli- 
able yet determined for these clusters. 

2. A significant number of relatively cool low gravity objects 
have been identified and are found on average to have large 
projected rotational velocity. The low gravities may then re- 
flect at least in part the significant centrifugal accelerations 
experienced by these stars. 

3. Observed Hertzsprung-Russell diagrams have been com- 
pared with theoretical predictions for both non-rotating and 
rotating models. The observed main sequence luminosities 



are in general larger than predicted. For NGC 6611, this 
may reflect non co-eval star formation. Other possible ex- 
planation include stellar rotation and the presence of a sig- 
nificant number of binary systems. 

4. Cluster ages have been deduced and some evidence for non- 
coeval star formation is considered especially for the young 
cluster NGC 6611. Indeed for none of the clusters does a 
single isochrone provide a good fit to the observed HR- 
diagram. 

5. Projected rotational velocities for targets in the older clus- 
ters, NGC 3293 and NGC 4755, have been found to be 
systematically larger than those for the field, confirming 
the conclusion of Strom et al. i2005j for h and x Persei. 
Contrary to the finding of Strom et al., this phenomenon 
appears to be present at all the B-type spectral types present 
in our sample. 

6. The distribution of projected rotational velocities are con- 
sistent with a Gaussian distribution of current rotational ve- 
locities. For the relatively unevolved targets in the older 
clusters, NGC 3293 and NGC 4755, the peak of the veloc- 
ity distribution would be 250km s^^with a full-width-half- 
maximum of approximately 180 km s^^. For NGC 6611, 
the corresponding values are 175km s^^and 165 km s^^. 

7. We find that the mean projected rotational velocity of stars 
which have strong winds (I5-6OM0) is lower than that of 
the lower mass stars. This is consistent with predictions of 
the rotating models which show the spin-down effect due 
angular momentum loss through stellar winds. 

8. For all three clusters we deduce present day mass functions 
with F-values in the range -1.5 to -1.8. 

Acknowledgements. We thank the referee, Steven Strom, for his valu- 
able and constructive comments on an earlier version of this paper. 
We are grateful the staff at both Paranal and La Silla for their in- 
valuable assistance with the observational programme. We also thank 
Jonanthan Smoker and Roheid Mokiem for discussions on rotational 
velocities and Georges Meynet for access to some of the low mass 
rotating models in advance of publication. CJE and JKL acknowl- 
edge financial support from the UK Particle Physics and Astronomy 
Reseai-ch Council (PPARC) under grants PPA/G/S/2001/00131 and 
PPA/G/O/2001/00473 and IH acknowledges financial support from 
the Northern Ireland Department of Education and Learning. AH, 
CAT and DJL thank the Spanish Ministerio de Educacion y Ciencia 
for support under project PNAYA 2004-0827 1-C02-01. SIS acknowl- 
edges support from the ESF in the form of a EURYl award. This pa- 
per is based in part on observations made with the Isaac Newton and 
William Herschel telescopes, operated on the island of La Palma by 
the Isaac Newton Group in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. 

References 

Abt, H. A., Levato, H., & Grosso, M. 2002, ApJ, 573, 359 
Afflerbach, A., Churchwell, E., & Werner, M.W. 1997, ApJ, 478, 190 
Balona, L.A. 1994, MNRAS, 268, 119 
Balona, L.A., & Koen, C, 1994, MNRAS, 267, 1071 
Baume G., Vasquez R.A., G. Carrai'o, & A. Feinstein, 2003, A&A, 
402, 549 

Belikov, A. N., Kharchenko, N. V., Piskunov, A. E., & Schilbach, E. 
1999, A&AS, 134, 525 



Dufton et al: Stellar parameters in NGC 3293, NGC 4755 and NGC 661 1 



17 



40 



25 



20 



1 5 



NGC 3293 
V„,= 300 km/s 










o 






'• 


oo 


• • 


o 


• 


• 




9 











Log (T.,) 



3 



40 



25 ° 

20 o 



15-" 8,,.--- 
12-' I 



NGC 4755 
V,„,= 300 km/s 



5-^ 



4,6 



4.4 4.2 
Log (T.„) 



4.0 



_ 60 o...-" 












o 



- 40 












-'• y 












25 












: 2o--' 


• 










! 15 ■ 


..-* 
•• • 












1 2 • 


• 










9 


o 








" NGC 661 1 

_ V,„,= 300 km/s 




t ° ... 
• 

• 

5 


o 

• o 

• 

. .-■ • 

o 

o 







Log (T.„) 



3 



Fig. 13. Hertzsprung-Russell diagrams for the three clusters together with rotating stellar evolutionary models from the Geneva 
group. The symbols are as in Fig|6l and all the models have an initial equatorial rotational velocity of 300km s^^. The tracks for 
the 25-6OM0models have been truncated before the bluewards evolution towards the WR phase, to aid clarity. 



18 



Dufton et al.: Stellar parameters in NGC 3293, NGC 4755 and NGC 661 1 



Bernacca, P.L., & Perinotto, M. 1974, A&A, 33, 443 

Blecha, A., North, P., Royer, R, & Simond, G. 2003, BLDR Software 

- Reference Manual, 1st. edition 

Brown, A. G. A., & Verschueren, W. 1997, A&A, 319, 81 1 
Bouret, J.-C., Lanz, T, Hillier, D.J., et al. 2003, ApJ, 595, 1182 
Chandrasekhar, S., Munch, G. 1950, ApJ, 111, 142 
Chini, R., & Wargau, W.F. 1990, A&A, 227, 213 
Crowther, PA., Hillier, D.J., Evans, C.J., et al. 2003, ApJ, 579, 774 
Crowther, PA., Lennon, D.J., & Walborn, N.R. 2006, A&A, 446, 279 
Daflon, S., & Cunha, K. 2004, ApJ, 617, 1115 
Dufton, PL., Smartt, S.J., & Hambly, N.C. 1999, A&AS, 139, 231 
Ehas, J. H., Frogel, J. A., & Humphreys, R. M. 1985, ApJS, 57, 91 
Evans, C.J., Crowther PA., Fullerton A. W., Hillier D. J., 2004, ApJ, 
610, 1021 

Evans, C.J., Smartt, S.J., Lee, J.K. & et al. 2005, A&A, 437, 467 
Evans, C.J., et al. 2006, A&A, submitted 

Dufton, P.L., Ryans, R.S.I., Trundle, C., Lennon, D.J., Hubeny, I., 

Lanz, T. & Allende Prieto, C. 2005, A&A, 
Feast, M. W. 1958, MNRAS, 118, 618 
Feast, M. W. 1963, MNRAS, 126, 111 

Freyhammer, L.M., Hensberge, H., Sterken, C., et al. 2005, A&A, 429, 
631 

Gies, D.R., & lambert, D.L. 1992, ApJ, 387, 673 
Girardi, L., Bressan, A., Bertelh, G., & Chiosi, C. 2000, A&AS, 141, 
371 

Gray D.F., 1992, The observation and analysis of stellar photospheres. 

Cambridge Univ. Press, 2nd ed., Cambridge. 
Gummersbach, C. A., Kaufer, A., Schaefer, D. R., Szeifert, T, & Wolf, 

B. 1998, A&A, 338, 881 
Heger, A., & Langer, N. 2000, ApJ, 544, 1016 
Herrero, A., Puis, J., Najarro, F 2002, A&A, 396, 949 
Hillenbrand, L.A., Massey, P, Strom, S.E., & Merrill, K.M. 1993, AJ 

106, 1906 

Hillier, D.J.,Lanz, T., & Heap, S.R., et al. 2003, ApJ, 588, 1039 
Hjorth, J., Sollerman, J., MoUer, P., et al. Nature, 423, 847 
Huang, W. Gies, D.R, 2005, 2005, ApJ, submitted, astro-ph/05 10450 
Howaith, 1. D., Siebert, K. W., Hussain, G. A. J., & Prinja, R. K. 1997, 

MNRAS, 284, 265 
Howarth, I.D., & Smith, K.C., 2001, MNRAS, 327, 353 
Hubeny, I., 1988, Computer Physics Comm., 52, 103 
Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 

Hubeny, I., Heap, S.R., & Lanz, T., 1998, in ASP Conf. Series 

- Boulder-Munich: Properties of Hot, Luminous Stars, ed. 
I.D. Howarth, 131, 108 

Hunter, 1., Dufton, PL., Smartt, S.J. et al, 2006, A&A, in preparation 
Kaufer, A., Stahl, O., Tubbesing, S., & et al. 1999, The Messanger, 
95, 8 

Kaufer, A., Stahl, O., Prinja, R. K., & Witherick, D. 2006, A&A, 447, 
325 

Keller, S.C. 2004, Publ. Astron. Soc. Austraha, 21, 310 

Kudritzki, R.-R, & Puis, J. 2000, ARA&A, 38, 613 

Kurucz, R. L. 2005, Memorie della Societa Astronomica Italiana 

Supplement, 8, 14 
Lynga, G. 1987, 10 European Regional Astronomy Meeting of 

the lAU, Prague, Czechoslovakia, Proceedings. Volume 4 

(A89-32590 13-90). Ondrejov, Czechoslovakia, Czechoslovak 

Academy of Sciences 121 
Maeder A., & Meynet G., 1991, A&AS, 89, 451 
Maeder A., & Meynet G., 2000, A&A, 361, 159 
Maeder A., & Meynet G., 2001, A&A, 373, 555 
Martayan C, Rremat, Y, Hubert, A. M., Ploquet, M., Floquet M., 

Zorec, J., Neiner, C, 2006, A&A, astro-ph/0601240 
Massey, P., Johnson, K. E., & De gioia-Eastwood, K. 1995, ApJ, 454, 

151 



McErlean, N. D., Lennon, D. J., & Dufton, R L. 1999, A&A, 349, 553 
Massey, R, Parker, J. W., & Garmany, C. D. 1989, AJ, 98, 1305 
Meynet, G., Maeder A., Schaller, G., Schaerer, D., & Charbonnel, G. 

1994, A&AS, 103, 97 
Meynet, G., & Maeder, A. 2000, A&A, 361, 101 
Meynet, G., & Maeder, A. 2003, A&A, 404, 975 
Mokiem, M.R., de Koter, A., Evans, C.J., & et al. 2006, A&A submit- 
ted 

Penny, L. R. 1996, ApJ, 463, 737 

Pasquini, L., Avila, G., Blecha, A., & et al. 2002, The Messenger, 1 10, 
1 

Puis, J., Urbaneja, M.A., Venero, R., & et al. 2005, A&A, accepted 
Rix, S., Pettini, M., Leitherer, C, et al. 2004, ApJ, 615, 98 
Rolleston, W. R. J., Smartt, S. J., Dufton, R L., & Ryans, R. S. I. 2000, 
A&A, 363, 537 

Ryans, R.S.I., Dufton, PL., Rolleston, W.J.R., Lennon, D.J., Keenan, 

FR, et al. 2002, MNRAS, 336, 577 
Ryans, R.S.I., Dufton, PL., Mooney, C.J., Rolleston, W.J.R., Keenan, 

FR, et al. 2003, A&A, 401, 1119 
Sagar, R., & Camion, R.D. 1995, A&AS, 111, 75 
Sanner, J., Bruzendorf, J., Will, J.-M., & Geffert, M. 2001, A&A, 369, 

511 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1 
Schaller G., Schaerer D., Meynet G., & Maeder A. 1992, A&AS 96, 
269 

Shobbrook, R.R. 1983, MNRAS, 205, 1215 

Slesnick, C.L., Hillenbrand, L.A., & Massey, R 2002, ApJ, 576, 880 
Smartt, S.J., Maund, J.R., Hendry, M.A., et al. 2004, Science, 303, 499 
Strom, S.E., Wolff, S.C, & Dror, D.H.A., 2005, ApJ, 129, 809 
Townsend, R. H. D. 1997, MNRAS, 284, 839 

Trundle, C, Lennon, D.J., Puis, J., & Dufton, PL. 2004, A&A, 417, 
217 

Trundle, C, & Lennon, D.J. 2005, A&A, in press 
Vacca, W.D., Garmany, CD., & Shull, M.J. 1996, ApJ, 460, 914 
Valquez, G.A., Leitherer, C, Heckmann, T.M., et al., 2004, ApJ, 600, 
162 

Vink, J.S., de Koter. A., Earners, H.J.G.L.M. 2001, A&A, 269,574 
Walborn, N.R., Morrell, N.I., Howarth, I.D., et al., 2004, ApJ, 608, 
1028 

Wolff, S.C, Edwards, S., & Preston, G.W. 1982, ApJ, 252, 322 
Wolff, S.C, Strom, S.E., Drod, D., Lanz, L, & Venn, K., 2006, ApJ in 
press 



Dufton et al.: Stellar parameters in NGC 3293, NGC 4755 and NGC 661 1 



19 



Fig. 3. Observed and rotationally broadened theoretical he- 
lium line profiles. 4755-033 (upper panel) and 4755-030 (lower 
panel) have estimated values of v sini of 75 and 330 km s~^, 
respectively. 
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Fig. 5. The procedure to correct the atmospheric parameter es- 
timates for stellar rotation is illustrated for a value of v sin i 
of 300 km/s. Plus sign {+) mark the atmospheric parameters 
deduced when the effects of rotation are included. Diamonds 
symbols (o) connected by a line, represent the original atmo- 
spheric parameters, which are designated by the adjacent pairs 
of numbers (note T^g is in x 1000 K). The correction factors 
are given in parenthesis. 
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Fig. 4. Initial estimates of atmospheric parameters for targets in 
NGC 3293, deduced from the hydrogen and helium equivalent 
widths. 
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Fig. 7. Present day mass functions for the three clusters. 



Fig. 8. Histograms of projected rotational velocities for each of 
the three cluster. 
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Fig. 11. The histogram of the observed projected rotational ve- 
locities for targets in the clusters NGC 3293 and NGC 4755 
is compared with models assuming an intrinsic rotational ve- 
locity distribution of a delta-function (a) and a uniform dis- 
tribution (b). The histogram data includes only stars in the 
mass range 3-15Mo i.e. the evolved objects have been re- 
moved. In panel (a) all stars have been assumed to rotate at 
velocities of 350km s^^(solid curve) and 250km s"^(dashed 
curve). In panel (b) continuous distributions of velocities be- 
tween zero and a maximum of 350km s~^(solid curve) and 
250km s^^(dashed curve) are assumed. It is clear that neither 
rotational velocity distribution provides convincing agreement 
with the observed velocity histograms. 
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Fig. 12. Upper Panel: The histogram of observed projected ro- 
tational velocities for targets in the clusters NGC 3293 and 
NGC 4755 is compared with the model assuming a Gaussian 
distribution of rotational velocities. The Gaussian has a maxi- 
mum at a rotational velocity 250 km s^^and a (T=110km s^^. 
The histogram data includes only stars in the mass range 3- 
15Mq i.e. the evolved objects have been removed. Lower 
Panel: The NGC 66 11 stars, with a Gaussian distribution with 
a maximum at 175km s^^and a o'=100km s^^ 
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Table 7. Physical parameters for NGC 3293: V, B ~V and spectral types are from Paper I; Atmopheric parameters (Tcff , log g) and projected 
rotational velocities (v sin i) were estimated using the methodologies discussed in Sect.|3| The specific method used to detemine the atmospheric 
parameters are designated as follows — He: from neutral helium and hydrogen equivalent widths; He+: from ionized helium, Silll/SilV 
ionization equilibrium and hydrogen line profiles; A: Tcff adopted; I: Teg interpolated. The radial velocities (Vr ) are taken from Paper I, whilst 
the distances from the cluster centre, r', are in arcminutes. 
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Table 7. Physical parameters for NGC 3293: V, B ~V and spectral types are from Paper I; Atmopheric parameters (Tcff , log g) and projected 
rotational velocities (v sin i) were estimated using the methodologies discussed in Sect.|3| The specific method used to detemine the atmospheric 
parameters are designated as follows — He: from neutral helium and hydrogen equivalent widths; He+: from ionized helium, Silll/SilV 
ionization equilibrium and hydrogen line profiles; A: Teff adopted; I: Teg interpolated. The radial velocities (Vr ) are taken from Paper I, whilst 
the distances from the cluster centre, r', are in arcminutes. 
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-15 


6.53 


3293-115 


13.11 


0.14 


B8 III-V 


12700 


4.10 


He 


170 


0.24 


2.16 


3.3 


-13 


1.73 


3293-116 


13.12 


0.11 


B6-7 V 


13000 


3.90 


He 


260 


0.21 


2.14 


3.4 


-14 


1.81 


3293-117 


13.15 


0.14 


B8 III-V 


13200 


3.95 


He 


34 


0.24 


2.19 


3.4 


-8 


3.67 


3293-118 


13.17 


0.06 


B8 III-V 


12100 


3.95 


He 


205 


0.15 


1.98 


3.1 


-10 


10.06 


3293-120 


13.20 


0.33 


B5 V 


14600 


4.10 


He 


245 


0.45 


2.52 


4.0 


-12: 


10.11 


3293-123 


13.26 


0.20 


B8 III 


12400 


3.85 


He 


275 


0.29 


2.14 


3.2 


-27 


1.09 


3293-124 


13.28 


0.21 


B8 III 


13200 


3.95 


He 


235 


0.31 


2.22 


3.5 


-9 


8.29 


3293-125 


13.29 


0.10 


B8 III-V 


12300 


3.95 


He 


210 


0.19 


2.00 


3.1 


-14: 


0.72 
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Table 8. Physical parameters for NGC 4755. The columns and sources of information are as for TableQ 



ID 


V 


U — V 


Sp. Type 


Jeff 


log 9 


Method 


V^sin? 


LLiyLi — V ) 


Inn T 

log Li 


Mass 


f r 


T 




mag. 






K 


-2 

cm s 




km s 








km s 




47'ss-nni 

"-r 1 O'J \J\J 1 


05.77 


0.32 


B9 la 


13500 


2.00 


He 


38 


0.44 


5.04 


20 


— 16 


2.85 


4755-002 


05.98 


0.22 


B3 la 


16300 


2.25 


Si 


70 


0.37 


5.18 


23 


— 10 


1.53 


4755-003 


06 80 


0.24 


B2 III 


24400 


3.40 


He 


38 


0.45 


5.35 


30 


— 21 


3.01 


4755-004 


06.92 


0.20 


B 1 5 lb 




2 75 


Si 


75 


0.38 


5.02 


20 


— 19 


0.32 


4755-005 


07.96 


0.20 


B2 III 


1 6900 


2.80 


He 


110 


0.36 


4.41 


12 


—53 


2.05 


4755-006 


08.37 


0.14 


Bl III 


21500 


3.15 




100 


0.33 


4.45 


13 


— 20 


2.68 


4755-007 


08.58 


0.1 1 


Bl V 


25000 


3.70 




170 


0.33 


4.52 


14.5 


— 12 


1.88 


4755-008 


08.61 


0.14 


BO. 5 V 


27000 


3.65 




170 


0.37 


4.61 


16.1 


—2 


0.48 




09.01 


0.1 1 


Bl V 


25000 


3.90 




115 


0.33 


4.34 


13 


— 15 


4.00 




09.38 


0.16 


Bl V 


25000 


4. 10 




180 


0.38 


4.26 


12.4 





1.71 




09.58 


0.17 


Bl 5 V 


22500 


4 10 


I 


22 


0.37 


4.06 


10 4 


—24 


1.12 


4755-012 


09.60 


0.16 


B1.5 V 


22500 


4.00 


J 


140 


0.36 


4.04 


10.3 


— 55 


1.81 


4755-013 


09.68 


0.15 


B1.5 V 


22500 


3.90 


I 


190 


0.35 


4.00 


10.0 


_3 


2.13 


'-r 1 JO -\J 1 J 


09 74 


0.16 


Bl V 


25300 


4.00 


Si 


48 


0.38 


4.13 


1 1.6 


— 18 


2.76 


4755-016 


09.76 


0.12 


B1.5 V 


22500 


3.80 


I 


195 


0.32 


3.93 


9.7 


— 8 


1.00 


47SS-01 7 


09.90 


0.18 


Bl 5 V 


22500 


4 10 


I 


75 


0.38 


3.95 


9.8 


— 15 


2.94 


4755-01 Q 


10.04 


0.17 


B1.5 V 


22500 


4.10 


I 


80 


0.37 


3.88 


9.4 


—20 


5.57 




10.05 


0.13 


B2 V 


19500 


3.80 


He 


38 


0.30 


3.65 


7 7 


—27 


6.68 


4755-021 


10.29 


0.20 


B1.5 Vn 


22500 


3.80 


I 


300 


0.40 


3.82 


9.2 


— 11 


3.26 


/ JJ~\JZ^Z^ 


10.89 


0.12 


B2 5 Vn 


18800 


3.70 


He 


300 


0.28 


3.26 


6.3 


— 18 


1.83 




10.99 


0.18 


B2 V 


19800 


3.95 


He 


120 


0.35 


3.35 


6.8 


—26 


1.85 


4755-025 


11.02 


0.14 


B2.5 V 


1 9600 


3.95 


He 


235 


0.31 


3.28 


6.6 


—20 


4.75 


/ JJ'\jA\j 


1 1. 19 


0.24 


B2 5 V 


1 8600 


3.95 


He 


205 


0.40 


3.27 


6.3 


— 15 


0.97 


47SS-097 


1121 


17 


B2 5 Vn 


1 7700 


3.80 


He 


210 


0.32 


3. 1 1 


5.7 


— 18 


1.05 


4755-029 


1 1.23 


0.19 


B2 V 


20200 


4.00 


He 


190 


0.37 


3.30 


6.8 


— 19 


3.02 


4755-030 

•-r / JJ~\JJ\J 


1 1.31 


0.21 


B2.5 Vn 


1 7500 


3.70 


He 


330 


0.36 


3.11 


5.6 


— 35 


3.66 


47SS-031 


1 1.34 


0.24 


B2 5 V 


1 8000 


3.95 


He 


150 


0.40 


3 17 


5.9 


— 20 


7.13 


47SS-039 


1 1.35 


0.19 


B2 5 V 


1 9000 


4.00 


He 


195 


0.36 


3. 17 


6.1 


— 13 


3.68 


4755-033 


1 1.35 


0.24 


B3 V 


17600 


4.00 


He 


75 


0.39 


3.14 


5.7 


— 19 


1.76 


4755-034 


11. 39 


0.23 


B3 V 


15200 


3.85 


He 


265 


0.36 


2.92 


4.8 


— 19 


2.91 


4755-035 


1 1.41 


0.23 


B5 V 


15100 


3.95 


He 


190 


0.35 


2.91 


4.7 


—31 


0.11 


4755-036 


11. 45 


0.16 


B3 Vn 


1 7500 


3.75 


He 


330 


0.31 


2.99 


5.4 


—29 


2.40 


4755-037 


1 1.48 


0.13 


B2.5 V 


18800 


4.05 


He 


195 


0.29 


3.03 


5.8 


— 19 


3.74 


47SS-038 


1 1 54 


0.16 


B3 Ye 


15400 


3 60 


He 


115 


0.29 


2.79 


4.6 


— 19 


2.72 


47SS-03Q 


1 1.54 


0.16 


B2 5 V 


18000 


4.00 


He 


235 


0.32 


2.99 


5.5 


—23 


1.24 


47SS_040 


1 1.54 


0.32 


B2.5 V 


19200 


4. 10 


He 


65 


0.49 


3.27 


6.4 


—20 


10.04 


47SS-041 

*-rlJJ -\jM- 1 


1 1.58 


0.21 


B3 V 


17100 


4.00 


He 


165 


0.36 


2.97 


5.3 


— 16 


0.96 


47SS-049 


1 1.58 


0.29 


B3 Vn 


17300 


3.80 


He 


345 


0.44 


3.08 


5.5 


— 14 


3.11 


4755-043 


1 1 59 


0.22 


B3 V 


1 6700 


4 10 


He 


125 


0.36 


2.95 


5.1 


— 29 


1.48 


4755-044 


11. 59 


0.28 


B5 V 


14300 


3.60 


He 


265 


0.40 


2.83 


4.4 


— 12 


8.80 


47SS.045 


1 1.61 


0.22 


B3 V 


17300 


3.95 


He 


245 


0.37 


2.99 


5.3 


— 17 


1.25 


47SS-04fi 


1 1.62 


0.20 


B3 Vn 


16300 


3.65 


He 


325 


0.34 


2.88 


4.9 


—29: 


3.80 


4755-047 


11.62 


0.23 


B3 V 


17100 


3.90 


He 


280 


0.38 


2.98 


5.3 


— 20 


1.97 


47'SS.n48 


1 1.63 


0.19 


B3 V 


1 9000 


4.15 


He 


55 


0.36 


3.06 


5.9 


— 20 


4.42 


47^5-049 


11. 65 


0.34 


B5 V 


14700 


4.00 


He 


140 


0.46 


2.91 


4.7 


— 36: 


8.46 


4755-051 

*-!■ / JJ~\JJ 1 


1 1.75 


0.23 


B3 Vn 


16200 


3.75 


He 


320 


0.37 


2.86 


4.9 


— 28 


2.25 


4755-052 


1 1.77 


0.17 


B3 V 


16600 


4.00 


He 


165 


0.31 


2.81 


4.9 


— 20 


4.71 


4755-053 


1 1.78 


0.18 


B3 Vn 


16000 


3.80 


He 


305 


0.31 


2 77 


4.7 


— 22 


2.96 


4755-054 


1 1 78 


0.19 


B3 V 


15300 


3.45 


He 


1 10 


0.33 


2.75 


4.7 


— 20 


3.98 


4755-056 


1 1.86 


0.29 


B3 Vn 


15600 


3.75 


He 


300 


0.42 


2.84 


4.7 


— 22 


1.25 


47^5-057 


11.94 


0.36 


B6-7 Ill-Ve 


13200 


3.65 


He 


235 


0.46 


2.69 


4.0 


—44: 


8.45 


/ JO \jJO 


12.07 


0.26 


B3 Vn 


15200 


3 75 


He 


320 


0.39 


2.69 


4.4 


—28 


0.98 


47SS-0(iO 

/ JO-\j\J\J 


12. 12 


0.19 


B5 V 


15400 


4.05 


He 


145 


0.32 


2.60 


4.3 


— 19 


3 77 


47SS-0fil 


12. 15 


0.25 


B5 V 


14500 


3.30 


He 


5 


0.39 


2.61 


4.3 


—20 


3.73 


4755-062 


12.18 


0.28 


B6-7 III-V 


14000 


4.00 


He 


200 


0.39 


2.57 


4.0 


— 14 


3.13 


4755-063 


12. 18 


0.33 


B8 V 


12400 


4.20 


He 


235 


0.42 


2.49 


3.6 


— 32 


5.72 


4755-064 


12.22 


0.27 


B6-7 V 


13500 


4.30 


He 


160 


0.38 


2.50 


3.8 


— 23 


1.05 


4755-065 


12.25 


0.18 


B3 V 


15200 


3.70 


He 


105 


0.31 


2.52 


4.2 


—22 


5.82 


/ JO \1\J\J 


12.25 


0.27 


B3 V 


15500 


4 10 


He 


165 


0.40 


2.65 


4.4 


—24 


1121 


4755-067 


12.28 


0.24 


B5 V 


14000 


4.05 


He 


290 


0.35 


2.48 


3.9 


— 23 


2.73 


4755-068 


12.31 


0.27 


B3 Vn 


14800 


3.95 


He 


270 


0.39 


2.57 


4.1 


— 20 


0.22 


4755-06Q 
*-!■ / JO \j\jy 


12.31 


0.28 


B5 V 


15000 


4.10 


He 


185 


0.40 


2.60 


4.2 


— 17 


2.16 


^ 1 JO -V / \j 


12.32 


0.23 


B3 V 


15200 


3.60 


He 


175 


0.36 


2.55 


4.2 


— 24 


3.75 


47SS-079 

*-r 1 JO \J 1 ^ 


12.43 


0.24 


B5 V 


14200 


4 10 


He 


100 


0.35 


2.44 


3.9 


— 19 


1.59 


4755-074 


12.47 


0.29 


B5 III-V 


13800 


3.50 


He 


120 


0.40 


2.45 


3.8 


— 36: 


2.58 


4755-075 


12 52 


24 


B5 III-V 


14100 


3 55 


He 


1 15 


35 


2 39 


3 8 


— 22 


1 71 


4755-076 


12.58 


0.36 


B5 V 


13800 


3.95 


He 


190 


0.47 


2.49 


3.8 


— 18 


0.43 


4755-077 


12.59 


0.32 


B6-7 III 


12800 


3.25 


He 


170 


0.44 


2.38 


3.8 


—27 


2.26 


4755-078 

^ 1 JO \J 1 o 


12.60 


0.32 


B5 Vn 


14200 


3.90 


He 


315 


0.43 


2.47 


3.9 


—23 


9.62 


4755-07Q 

"-T 1 JO \j 1 y 


12.62 


0.28 


B5 V 


13700 


3.85 


He 


250 


0.39 


2.37 


3.7 


— 20 


2.41 


4755-080 


12.62 


0.35 


B8 III-V 


12200 


3.95 


He 


220 


0.44 


2.32 


3.4 


— 19 


3.86 


4755-08 1 


12.63 


0.28 


B6-7 Illn 


13200 


3.65 


He 


285 


0.38 


2.32 


3.5 


— 18 


1.67 


4755-082 


12.65 


0.30 


B8 III-V 


12000 


3.80 


He 


245 


0.39 


2.22 


3.2 


-25 


2.93 


4755-087 


12.81 


0.25 


B8 III 


13500 


4.10 


He 


280 


0.36 


2.24 


3.5 


-19 


1.18 


4755-088 


12.83 


0.25 


B8 Illn 


12900 


3.80 


He 


260 


0.35 


2.17 


3.4 


-24 


6.34 


4755-089 


12.84 


0.36 


B8 III 


12600 


3.35 


He 


160 


0.48 


2.31 


3.6 


-19 


6.94 


4755-090 


12.88 


0.27 


B8 Illn 


12600 


3.65 


He 


290 


0.36 


2.15 


3.3 


-16 


2.39 


4755-091 


12.90 


0.37 


B8 III 


14100 


4.40 


He 


50 


0.48 


2.40 


3.8 


-20 


4.11 


4755-093 


12.91 


0.28 


B8 III-V 


13200 


4.35 


He 


190 


0.38 


2.21 


3.4 


-29 


1.09 


4755-094 


12.93 


0.31 


B8 III 


13400 


3.80 


He 


90 


0.41 


2.25 


3.5 


-30 


4.52 
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Table 8. Physical parameters for NGC 4755. The columns and sources of information are as for TableQ 
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ID 


V 


U — V 


Sp. Type 


Jeff 


log 9 


Method 




LLiyLi — V ) 


Inn T 

log Li 


Mass 




T 




mag. 






K 


-2 

cm s 




, - 1 
km ,s 








1 - 1 

km s 




4755-095 


12.96 


0.20 


B8II1 


12700 


4.00 


He 


225 


0.30 


2.04 


3.2 


-22 


5.51 


4755-096 


12.96 


0.29 


B8 111 


13300 


4.25 


He 


180 


0.39 


2.21 


3.5 


-13 


1.48 


4755-098 


13.00 


0.25 


B8 Illn 


12700 


3.95 


He 


280 


0.35 


2.09 


3.3 


-20 


3.93 


4755-099 


13.01 


0.41 


B8 III 


12600 


3.75 


He 


250 


0.50 


2.27 


3.4 


-20 


9.25 


4755-100 


13.09 


0.26 


B8 III 


13300 


4.40 


He 


220 


0.36 


2.12 


3.4 


-15 


1.09 


4755-106 


13.16 


0.32 


B8 III 


12200 


3.80 


He 


290 


0.41 


2.06 


3.1 


-18 


9.89 


4755-107 


13.18 


0.29 


B8 Ill-V 


12600 


4.10 


He 


165 


0.38 


2.06 


3.2 


-20 


2.43 


4755-108 


13.22 


0.31 


B8 Ill-V 


12400 


4.05 


He 


165 


0.40 


2.05 


3.2 


-22 


3.70 
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Fig. A.l. VLT-FLAMES targets (open circles) and FEROS targets (open squares) in the NGC 661 1 field. 

Fig. A.2. VLT-FLAMES targets (open circles) and FEROS targets (open squares) in the inner 7' of the NGC 3293 field. 

Fig.A.3. VLT-FLAMES targets (open circles) and FEROS targets (open squares) in the NGC 3293 field. Most of the objects 
observed in the inner region are shown separately in Figure [A2l 

Fig. A.4. VLT-FLAMES targets (open circles) and FEROS targets (open squares) in the inner 7' of the NGC 4755 field. 
Table 9. Physical parameters for NGC 6611. The columns and sources of information are as for Table|7| 



ID 


V 


B -V 


Sp. Type 


T tt 


log 9 


Method 


Vsinz 


E{B - V) 


log L 


Ma,ss 




r 




mag. 






K 














km s^^ 




6611-OOlt 


08.18 


0.34 


BO III 


32000 


3.60 


He-^ 


142 


0.60 


5.39 


40 


86 


13.17 


6611-002 


08.18 


0.43 


04 Ill((f + )) 


41500 


3.90 


He-^ 


102 


0.71 


5.86 


63 


14 


1.29 


6611-003 


08.73 


0.45 


06-7 V((f)) 


40000 


3.90 


He-^ 


87 


0.73 


5.63 


50 


17 


1.08 


661 1-004 


08.90 


0.04 


08.5 V 


37000 


4.00 


He-^ 


76 


0.32 


4.85 


24.0 


23 


4.18 


661 1-005 


09.13 


0.44 


08 III 


36000 


3.80 


He-^ 


90 


0.71 


5.32 


25 


5 


29.49 


661 1-006 


09.39 


0.24 


09.7 IIIp 


34000 


4.10 


He-^ 


43 


0.51 


4.84 


21.9 


-1 


4.04 


6611-008 


09.46 


0.82 


07 11(f) 


36000 


3.50 


He-^ 


100 


1.09 


5.75 


60 


15 


1.79 


6611-010 


09.75 


0.49 


Bl: e 






He 


225 








4 


12.22 


6611-011 


09.85 


0.58 


09 V 


36000 


4.20 


He-l- 


66 


0.85 


5.24 


29.5 


17 


1.54 


6611-012 


09.85 


0.48 


B0.5 V 


27000 


4.00 


Si 


75 


0.71 


4.70 


17.0 


14 


3.75 


6611-015 


10.12 


0.43 


09.5 Vn 


32500 


3.90 


He-^ 


410 


0.69 


4.77 


20.3 


23: 


0.72 


6611-017 


10.37 


0.57 


09 V 


36000 


3.95 


He-^ 


95 


0.84 


5.02 


25.6 


19 


2.54 


6611-018 


10.56 


0.35 


B8 III 


12900 


3.90 


He 


90 


0.45 


3.28 




16 


9.34 


6611-019t 


10.68 


0.36 


B1.5 V 


22500 


3.85 


I 


235 


0.56 


3.96 


9.8 


-47 




6611-020 


10.69 


0.40 


BO. 5 Vn 


27000 


4.10 


A 


220 


0.63 


4.24 


13.0 


20 


6.17 


6611-021 


10.80 


0.47 


Bl V 


26300 


4.25 


Si 


32 


0.70 


4.27 


12.9 


-18 


0.30 


6611-022 


10.98 


0.82 


Herbig Be 






He 


375 










0.46 


6611-023 


10.99 


0.41 


B3 V 


14600 


3.35 


He 


165 


0.55 


3.38 




-19: 


7.37 


6611-025 


11.20 


0.59 


Bl V 


25000 


4.10 


A 


95 


0.81 


4.23 


12.2 


8 


0.73 


6611-027 


11.26 


0.45 


Bl V 


25000 


4.10 


A 


255 


0.67 


4.00 


10.9 





2.85 


6611-029 


11.29 


1.05 


08.5 V 


36000 


3.85 


He-^ 


135 


1.32 


5.37 


35 


-6 


4.53 


6611-032 


11.48 


0.36 


B1.5 V 


22500 


4.15 


I 


70 


0.56 


3.64 


8.5 


4 


7.46 


6611-033 


11.47 


0.56 


Bl V 


25600 


4.00 


Si 


30 


0.78 


4.11 


11.7 


8 


1.07 


6611-034 


11.47 


0.55 


B8III 


13600 


4.05 


He 


60 


0.66 


3.28 




-18 


7.76 


6611-035 


11.51 


0.75 


B0.5 V 


27000 


4.10 
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Appendix A: Finding chiarts for Gaiactic ciusters 



With the benefi t of hindsight on the part of the authors, we now also include finding charts for the targets in the three Galacti c clus ter s repo ited in Paper I. The NGC 66 1 1 field is sho wn 
in Figure cores of NGC 3293 and NGC 4755 are fairly dense, so we show the central l' of each cluster in Figure s |A.2| and |A!4| with the full fields shown in Figures |A.3| and 
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Fig.A.5. VLT-FLAMES targets (open circles) and FEROS targets (open squares) in the NGC 4755 field. Most of the objects 
observed in the inner region are shown separately in Figure lA!4l 
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